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Abstract. We prove that the critical Maxwell-Klein Gordon equation on M 4+1 is globally 
well-posed for smooth initial data which are small in the energy. This reduces the problem 
of global regularity for large, smooth initial data to precluding concentration of energy. 



1. Introduction 

Let M 4+1 be the five dimensional Minkowski space equipped with the standard Lorentzian 
metric g = diag(l, — 1, — 1, — 1, — 1). Denote by : M 4+1 — > C a scalar function, and by 
A a : M 4+1 — > R, a = 1 . . . , 4, a real valued connection form, interpreted as taking values in 
isu(l). Introducing the curvature tensor 

F a p '■— d a Ap — dpA a 

as well as the covariant derivative 

D a (j) := (d a + iA a )(f) 

the Maxwell- Klein Gordon system are the Euler-Lagrange equations associated with the 
formal Lagrangian action functional 

C(A a , 0) := ~ / (\D a( £D^4> + -F a pF aP ) dxdt; 
Z J R 4+i Z 4 

here we are using the standard convention for raising indices. Introducing the covariant wave 
operator 

□ A := D a D a 

we can write the Maxwell-Klein Gordon system in the following form 

d f) F aP = -J a :=$l(<f>D^), 
n A( f) = 

One key feature of this system is the underlying Gauge invariance, which is manifested by 
the fact that if (A a , (p) is a solution, then so is (A a — d a \, e zx (f>). This allows us to impose an 
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additional Gauge condition, and we shall henceforth impose the Coulomb Gauge condition 
which requires 

4 

(2) £a^ = ° 

i=i 

The MKG-CG system can be written explicitly in the following form 

(3a) UAi = ViJ x 

(3b) n A </> = 

(3c) AA = J 

(3d) Ad t A = V l J h 

where the operator V denotes the Leray projection onto divergence free vector fields, 

V = I-VA _1 V 

The last equation (13d|) is a consequence of (I3"cl) due to the divergence free condition on the 
moments d a J a = 0, which in turn follows from (13bl) . However, it plays a role in the sequel 
so it is added here for convenience. 

Here it is assumed that the data Aj(0, •) satisfy the vanishing divergence relation (j2J). We 
remark that given an arbitrary finite energy data set for the MKG problem, one can find a 
gauge equivalent data set of comparable size which satisfies the Coulomb gauge condition. 
This argument only involves solving linear elliptic pde's, and is omitted. The key question 
now is to decide whether a family of data 

(A, 0)[O] := (4,(0, •), d t A a (0, •), 0(0, •), d t cf>(0, •)) 

which satisfy the compatibility conditions required by the two equations for 4 above can 
be extended to a global-in-time solution for the Maxwell-Klein-Gordon system. The key for 
deciding this question is the criticality character of the system. Note that the energy 
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is preserved under the flow ([TJ , and in our 4 + 1-dimensional setting it is also invariant under 
the natural scaling 

(f)(t, x) -»■ X(j)(Xt, Ax), A a (t, x) -> \A(\t, Ax) 

This means the 4+1-MKG system is energy critical, and in recent years a general approach to 
the large data Cauchy problem associated with energy critical wave equations has emerged. 
The first key step in this approach consists in establishing an essentially optimal global well- 
posedness result for data which are small in the energy norm, which is usually the optimal 
small-data global well-posedness result achievable. In this paper, we set out to prove this 
for the 4 + 1-dimensional Maxwell-Klein-Gordon system. 

Theorem 1. a) Let (A,<f))[0] be a C°° Coulomb data set satisfying 

(4) £(40)<e* 

for a sufficiently small universal constant e* > 0. Then the system (CQ) admits a unique global 
smooth solution (A, <fi) on M 4+1 with these data. 

2 



b) In addition, the data to solution operator extends continuously^ on the set (111) to a map 
H^R 4 ) x L 2 (R 4 ) 3 (A0)[O] -> (A,(f>) G C(R; if^lR 4 )) fl C yl (IR, L 2 (M 4 )) 

We remark that the same result holds in all higher dimensions for small data in the scale 
invariant space if 2^ x H^~ 2 . This has already been known in dimensions n > 6, see (21 j . 
We have chosen to restrict our exposition to the more difficult case n = 4 in order to keep 
the notations simple, but our analysis easily carries over to dimension n — 5. On the other 
hand, we do not know whether a similar result holds in dimension n — 3. 

Before explaining some more details of our approach, we recall here earlier developments 
on this problem, and how our approach relates to these. Considering the case of general 
spatial dimension n and denoting the critical Sobolev exponent by s c = | — 1 (thus s c — 1 for 
n = 4, corresponding to the energy), a global regularity result for data which are smooth and 
small in H Sc was established in dimensions n > 6 in the work [21] which served as inspiration 
to the present work. We note in passing that a result analogous to [21] was established in 
[18] for the Yang-Mills system in dimensions n > 6, and the present work most likely also 
admits a corresponding analogue for the 4 + 1-dimensional Yang-Mills problem. The global 
regularity question for the physically relevant n = 3 case of the Yang-Mills problem had 
been established earlier in the groundbreaking work [6]. Observe that this problem is energy 
sub- critical. 

In the context of MKG, the result [21] had been preceded by a number of works which 
aimed at improving the local wellposedness of MKG in the n = 3 case, beginning with 
[9], followed by [5], and more recently [19]; the latter in particular established an essentially 
optimal local well-posedness result by exploiting a subtle cancellation feature of MKG, which 
also plays a role in the present work. We also mention the recent result [3J which establishes 
global regularity for energy sub-critical data in the n = 3 case, in the spirit of earlier work 
by Bourgain [3]. 

In the higher dimensional case n > 4, an essentially optimal local well-posedness result for 
a model problem closely related to MKG was obtained in [14J. This model problem docs not 
display the crucial cancellation feature of the precise MKG-system which enable us here to 
go all the way to the critical exponent and global regularity. We mention also that essentially 
optimal local well-posedness for the exact MKG-system was obtained in [17]. Finally, the 
recent work [24] established global regularity of equations of MKG-type with data small in 
a weighted but scaling invariant Besov type space, in the case n = 4. 

The present paper follows a similar strategy as [2T] : one observes that the spatial Gauge 
connection components A,, j = 1, 2, 3, which are governed by the first equation fl3aj) . may in 
fact be decomposed into a free wave part and an inhomogeneous term (the second term in 
the Duhamel formula for A) which in fact obeys a better ^-Besov type bound (while energy 
corresponds to I 2 ) 

yjfrec _|_ yjnonlin 

This is important for handling the key difficulty of the MKG-system , which is the equation 
for (p, i. e. the second equation fl3bl) . In fact, we shall verify that the contribution of the term 
A nonlm to the difficult magnetic interaction term 2iAjdj(f> in the low-high frequency interaction 
case can be suitably bounded when combined with the term 2iAodt<f), an observation coming 
from [19]. However, the contribution of the free term A iTCC to the magnetic interaction term 



here the continuity is locally in time 
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is nonperturbative and cannot be handled in this manner. Thus, following the example set 
in [21], we retain this term into the covariant wave operator^ A f ie e. More precisely, we shall 
define a suitable paradifferential wave operator D P A which incorporates the 'leading part' of 
□^frco while relegating the rest to the source terms on the right. 

The key novelty of this paper then is the development of a functional calculus, involving 
in particular X s ' b -type as well as atomic null-frame spaces developed in other contexts, for 
solutions of the general inhomogeneous 'covariant' wave equation 

* P A = f 

This refined functional calculus is necessary to control the nonlinear interaction terms, which 
become significantly more delicate in the critical dimension than in the setting studied in 
|21j . In particular, the Strichartz norms themselves appear far from sufficient to handle 
the present situation. The above covariant wave equation will be solved by means of a 
suitable approximate parametrix, and we show that this parametrix satisfies many of the 
same bounds as the usual free wave parametrix, in particular encompassing refined square 
sum type microlocalized Strichartz norms as well as null-frame spaces. We expect the calculus 
developed here to be of fundamental importance in other contexts, such as the regularity 
question of the critical Yang-Mills system and related problems from mathematical physics. 

2. Technical preliminaries 

Throughout the sequel we shall rely on Littlewood-Paley calculus, both in space as well 
as space-time. In particular, we constantly invoke the standard Littlewood-Paley localizers 
Pk, k £ Z, which are defined by 

Kf = x(§)f(0 

for functions / defined on IR 4 . Here x is a smooth bump function, supported on [~, 4], which 
satisfies the key condition ^ fceZ x(^) = 1 if £ > 0. To measure proximity of the space-time 
Fourier support to the light cone, we use the concept of modulation. Thus we introduce the 
multipliers Qj, j £ Z, via 

with the same \ as before, where "in this context denotes the space-time Fourier transform. 
We then refer to 2 J as the modulation of the function. On occasion we shall also use 
multipliers Si, which restrict the space-time frequency to size ~ 2 l . These multipliers allow 
us to introduce a variety of norms. In particular, for any p £ [1, oo), we set for any norm 

IMI«==(EM) i 

fcez 

We also have the following X s ' 6 -type norms, applied to functions localized to spatial fre- 
quency ~ 2 k : 

\\F\\ x; ,:= 2^(^[2^||g,.F|| L? J p )^£[l,oo), 



3 Hcre we set ,4f, rcc = 0. 
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with the obvious analogue 



|^:=2 sfe sup 2^1^.4 2 



t,x 



kez 

For more refined norms, we shall also have to use multipliers Pf, which localize the homo- 
geneous variable A to caps 00 C S 3 of diameter 2 l , by means of smooth cutoffs. In these 
situations, we shall assume that for each I a uniformly (in I) finitely overlapping covering of 
S 3 by caps u has been chosen with appropriate cutoffs subordinate to these caps. Similar 
comments apply to the multipliers P c v which localize to rectangular boxes and will be de- 
fined below. 

Given a norm || • ||s with corresponding space S, we denote by Sk the space of functions in 
S which are localized to frequency ~ 2 k . Furthermore, we denote by 

Sk,± 

the subspace of functions in Sk with Fourier support in the half-space r >< 0, with r the 
Fourier variable dual to t. 

3. Function Spaces 

There are three function spaces we work with: N, N*, and S. These are set up to that 
their dyadic subspaces Nk, N%, and S k satisfy the following relations: 

(5) N k = L\L 2 ) + Xl'~* , X°/> CS k C N* k , 

Then define: 

\\F\\ 2 N = EW P ^Wk- 
k 

We also define S k by 

\\ u \\s{ = \\ nu \\N k + ||Vm||l^l2 
On occasion we need to separate the two characteristic cones {r = ±|£|}. Thus we define 

N k>± , N k = N k ,+ n N k: _ 
ctl ctl ctt 1 at 

^* ± , m = N* k>+ + N* k> _ 

Our space Sk scales like L 2 free waves, and is defined by: 



|2 _ || M 11 2 1 II M 112 

OS 

'A- 



s k — 11 v \\s? tr + 11 v \\sr g 



where: 



(6) Sf = n i+ 3 Z2< 32^ + " 2 ^(L r ) , U\\ 2 = sup^||PrQ<fcWHko , 



The angular sector norms S k (l) are essentially the same as in the wave maps context, see 
e.g. [27], and defined shortly. Our space of solutions scales like free waves with H 1 data, 
with a high modulational gain as in [25]. Thus we set: 



(7) || = II V ^ Wsu + II D( P 



2 
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For later reference, we shall also use the norms 

?n-i 



Returning to S%(1), we define these as usual except that we need to add additional square 
information over smaller radially directed blocks Cfc'(Z') dimensions 2 k x (2 k +l ) 3 with ap- 
propriate dyadic gains. First define: 



NE 



SUP \\fu<t>\\L°?(L* ,) > 



where the norms are with respect to = t ± u ■ x and the transverse variable, while y ' u 
denotes spatial differentiation in the {£^) ± plane. Now set: 

(8) !I0||^ (O = H\\l r + 2- 2k u\\% E + 2- 3k Y,\\Q ± <P\\] wm 

± 



+ sup V (\\Pc k ,(i')<P\\ 2 s f r + 2" 2fc ||P c 



1 \\ne 

k i^l\> S7 (lf\ 



2- 2fe '-1 Pc AV) 4> \\b^) + 2" 3 ^ £ II Q ±P Cy 



± 



We remark that an important feature of these norms is that the time- like oriented L 2 (L°°) 
block norm gains dyadically from the length of Cfc'(Z') in the radial direction, while the null 
oriented L^(L^) norms gain from the size of Ck>(l') in the angular direction. We also remark 
that a useful feature of this setup is: 

(9) (EH P ^ P ^Hiw)" ^ 2 fc '2^||P fc 

where Cy are a finitely overlapping set of cubes of side length 2 k> . This follows by splitting 



0,4 1 
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Pk4> — Q<k'Pk4> + Q^k'Pk4> an d using the S% (\(k' — k)) norm for the first term and the X, 
norm and Bernstein's inequality for the second. 

Next we describe some auxiliary spaces of L l (L°°) which will be useful for decomposing 
the non-linearity. The first is for the hyperbolic part of the solution: 

|| <fi \\z h vp = / , II Pk<\> \ z hyv ' II^IIz^p = SU P / . 2'|| PfQk+2l4> IliifL 00 ) • 
, fe k KC 

k LO 

Note that as defined this space already scales like H 1 free waves. In addition, note the 
following useful embedding which is a direct consequence of Bernstein's inequality: 

(10) □- 1 £ 1 L 1 (L 2 ) C Z hyp . 

The second is for the elliptic part of the solution: 

ii^iu«" — ^2 wqj^Wlhl^) , \\A\\ Z di = 11 PkA \\ z m . 

j<k+C k 

We remark that the only purpose of this last norm is to handle sums over Qj in product 
estimates involving A in 
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Finally, the function spaces for Aq are much easier to describe as the Aq equation is elliptic: 

|2 _ II VT A l|2 , || A l|2 , lla A ||2 

4" 



Pollyi = || V^Aollioo^ + ||A || 2 2 . 3 + \\d t A \\ 2 2 . i 



We also have the derivative norms 

||Ao||y* := ||Vf-%||yi,iV>2 

4. Decomposing the non-linearity; Statement and use of the core 

multilinear estimates 



Recalling the definition of the currents J a = —5s((J)D a <j)) we write the MKG-CG system 
again here as: 

(Ha) UAi = ViJ x , 

(lib) u A< j> = , 

(11c) AA = J , 

This system will be solved iteratively using the following scheme. We initialize 

Af ] = Af ee , A 1} = 0, = frcc , 

where Af ee and free solve the flat wave equation with initial data Aj[0], respectively 0[O]. 
Given Aa , and their associated currents ji , we define the next iteration via the 
equations 

(12a) uA { ; n+l) = 7Vi m) , 

(12b) □ A(m) (m+1) = , 

(12c) AA m+1) = 4 m) , 

with the same initial data Aj[0], respectively 4>[0}. Assuming small energy for the initial data 
Ai[0] and 0[O] as in (J4J), we will prove that this Picard type iteration converges in the space 
S 1 . Our starting point is the linear bound 

(13) |Mj?V + ll* (1 V<Cbe. 
Then we will inductively establish the bound 

(14) ||4 m+1 ) - A^\, s , + ||0( m+1 > - <P^\\ S , + ||4 m+1) - 4 m) || yl < (Ce*)™ 

for a universal constant C > 2 Co- 

Assuming this holds, passing to the limit as n — > oo we obtain a Coulomb solution (A, <fi) 
to the MKG equation which satisfies the bound 

(15) IK° nlin H^ + U\W + Poll^ <e* 

The same argument proves uniqueness. If two solutions (A^,^ ^) and (A^\ (p^) have the 
same Cauchy data for A x , then the same A^ ee is used for both. Thus applying the same 
series of estimates in this section to the difference of the two exact solutions rather than two 
approximate solutions we obtain the bound 

(16) ||4 0) ~A X %^ + ||0 (O) -0 (1) ||si + ||4 0) -A^\\y, < ||0 (O) [O] -^[0]|| HlxL2 
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This bound proves both uniqueness and Lipschitz dependence of the solution with respect 
to 0[O]. The continuous dependence with respect to A x [0] is a more delicate issue and will 
be explained in section [5j 

The estimate fll4p (m) will follow from f )14p (< m). Summing up f)14p (< m) and (TT5|) . we 
can easily add to our induction hypothesis the bound 

(17) ||4 B) - + U {n) \\s^ + H4 n) Hi- < 2C e*, n<m 

provided that e* is small enough. 

A simple but very useful observation is that, since all iterated solve covariant wave 
equations, it follows that the associated moments are divergence free, D a A m) = 0. Then, 
differentiating the equation (I12cl) . we obtain 

(18) A^4 m+1) = d\j\ m) 

which will be used to estimate the high modulations of A . This is the reason why we have 
completely avoided using </>( m ) in (112bp . even though some of the terms in there will be 



treated perturbatively. 

A second fact to keep in mind is that while in terms of formulas this is a one step iteration, 
in terms of estimates this is really a two step iteration. Precisely, in order to obtain good 
bounds for <^>( m+1 ) we will need to reiterate "bad" portions of A a in terms of </>( m-1 ) (but not 
A*-" 1-1 )). All this is explained in detail below. 

To decompose the non-linearity, the following device will be useful. If Ai(D t>x , D t y ) is any 
bilinear operator we set: 

U k M{<p^) = QjPkM(Q <j - C <l>,Q< j -cil>) , 

j<k+C 

HtM(4>,i>) = QKj-cMiQjPkhQ^-crf) , 

j<k+C 

Now we describe the function spaces and multilinear estimates for each iterative piece of 
the non-linearity. In this Section we only state the main estimates, and reduce them to 
appropriate dyadic bounds. These will be proved later in the paper. 

4.1. Estimates for the A{. We split the spatial potentials into homogeneous and inhomo- 
geneous parts with respect to t = as follows, A\ m ^ = Af oc + ^ onlm >( m ) w here the second 
part solves the linear equation 

(19) □ A nonlin,(m+l) = ^(m) y^m) + j^M^M^ A™ nlin [0] = 

In order to establish the / 1 5' 1 bound in the first term of f[T4l it suffices to work directly with 
this equation. However, A^ also appears in the 0( m+1 ) equation f!12bl) . and in order to be 
able to treat its contribution perturbatively there we also need to control its Z hyp norm (see 
the subsection below devoted to (p). This works for the most part, but there is a part of 
requires a more delicate treatment, which is: 

(20) ^nonWH = _ g □- 1 ^(^ 1) V,^- 1) ) • 

A;<min{fci,fc2}— C 
8 



For the good part of J 4° onlm ' (m ) we w [\\ establish the difference bound 
(21) 

< (Ce*) n ~\ 3<n< m+1 

as well as the summed estimate 

(22) || A nonlin,(n) _ ^nonlin.W ^ < ^ 2 < 72 < m + 1 

On the other hand the bad part ^ J 4^ lonlm '( m ) considered later in combination with the 

similar part of Aq. 

As seen above, we need estimates not only for the equation ( fl9i) . but also for differences of 
two consecutive iterations. Thus we are led to consider a more general equation of the form 

(23) UBi = V0 l) V x (t) {2) + (3) (4) 4 (1) ) , 5™ on/m [0] = 0. 
and its corresponding bad part, 

(24) UBi = - u ~ X UkVM}VA { £) ■ 



fc<min{fci ,k,2}—C 



Then we will show: 



Proposition 4.1. One has the following iterative estimates for the functions B^ and JiBi 
defined above: 

(25) || £,11*51 < ll^llsill^llsi + IIA^II^II^Ikll^lUi, 

(26) \\B t -UB t \\ zhyP < || |M| (2) lis 1 • 

In combination with our induction hypothesis, the estimate (125]) leads to the / 1 5' 1 bound 
for the first term in (JHJ). Similarly, the estimate (126]) yields the Z hyp bounds for the good 
part of 4 nonlin ' (m) in ^ and 

Decomposition of the proof of estimates ( 125]) and ( 126]) . For estimate ( 125]) we begin with the 
high modulational bounds: 

(27) ||0 (1) V^ ll^-i, < H {1) \\s4<P {2) \\si , 

(28) ||WV 5) \\ elL2{A - h < ll^ll^ll^lkll^ll^. 



Next, for the cubic terms in A™ onhn both ( [25]) and ( 126]) follow once we can show: 

(29) || ||* L1(i2) < || ^ \\s4 || sl || ^ \y . 

Returning to the quadratic part, we employ the notation: 

Mj(0 (1) ,0 (2) ) = drfWdrfV - d^d4 m , 
which allows us to write: 

(30) V0Vv x <j> m ) = A^V%(0«,0( 2 )) . 
Then for the quadratic part of (125]) it suffices to show: 



(31) || V*A-W(0« (2) )||* JV < || Ml (2) |k 
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where M denotes any instance of the J\fij. Finally, making an analogous definition to (jS 
for each term in fl30l) we need to show: 



(32) 



(J-W)V S A-W(^ 



\( 2 )\ ii < 

; ) \\UZ h VP 



a 2 ) 



□ 



4.2. Estimates for Ao and OiAq. The estimates for temporal potentials are analogous to 
those above, but using instead the equations (I12cl) and (ITSj) . In expanded form these are 
written as 



AA 



(m+l) 



-9 



Ad t A { ™ +1) = -d l ^ {m) di^ m) ) - 9 4 (4 m) |0 (m) | 2 ) 

In a first approximation, from Aq"^ we isolate to output of bilinear high x high — > low 
interactions, namely 



-4 



(m+l),hh 




E 



fc<min{fci,A:2}— C 



To bring the analysis to the same level as in the case of the A/s, in analogy with (|20|) . we 
also define the component 7iA^ l+1 ^ of J \^ n+1 ^ hh j 



(33) 



(m+l) 



E 



( m )^ A m )> 



fc<min{fci ,k2]—C 



In addition to the estimates included in f ll4p . we also need the some similar bounds for 
j^(m+i),hh) as we n as bounds for the good differences for 3 < n < m + 1: 



1-4 



(n),hh 



A 



t 1)M \\ elLHHh < (CeT-i 



(34) 



||(A (n) _ A (n), hh) _ {A (n-l) _ < 

iK^4 n) - 4 n)M ) - (^4 n " 1} - 4 n_1),hh )iu«« < (Ce)"- 1 

Passing to differences we arrive at a system of equations with multilinear inhomogeneities 
of the form: 



(35) AB 

(36) Ad t B 
As above we write: 



0(D^0(2) +0 (3) (4) A (D ^ 

v,(0«v^ + 0( 3 V (4) A (1) ) • 



B, 



hh 



E 

k < m i n { k i , /c 2 } — C 



respectively 

(37) HB 
Then we will show: 



E 



fc<min{fci,A;2}— C 



10 



Proposition 4.2. One has the following estimates for Bo and dtBo defined above: 



(38) 


1 V^je-Bo | L°°(L 2 ) 


< 1 


0« 


S 1 


</> (2) II 


(39) 


II (B ,B% h ) \\ ilL2{A ^ 


< 1 

r^j 1 


0« 


Is 1 ! 


{2) II 


(40) 


II 9tB H^a^j 


< 1 


0« 


S 1 


(2) II 


(41) 


II -Bo — -Bo^ II^l 1 ^ 00 ) 


< 1 






(2) II 


(42) 


-Bp'" — H-B \\z m 


< 1 


0« 


U 1 ! 


4> {2) II 



s 1 + ii 4 1} 



s 1 + ||A W lis 



li^f/ 1 ) 

'l2(£T§)!I 
,0) 



(3) Us 1 

/A(3) II „, 



.4 



(1) II 3 I 

»iP(Hiy 



I SMI 



6^ 
Is 1 

AW 
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The bound for the last term in ( fl4l) follows from ( 1381) . ( 139]) and ( 1401) . The three bounds 
in ()34]) are consequences of (i55|) . ( 1411 and ( l42il . 

Outline of proof of estimates f p9|) -f l4Tl) . Estimate fl38|) follows from Sobolev's embedding, 
while the proofs of fl39|) and (140"]) either follow immediately from ( |2"T|) and fl28|) above, or in 
the case of Ao h from the estimates used to produce those bounds. Estimate ( I4TI) follows from 
similar considerations. In detail, the estimates ( 138]) - ( I4T]) will be proved in subsection 14.41 
On the other hand (T4"2"]) is more microlocal in nature and follows from calculations similar to 
those in the proof of fl3"2~]) listed above. It will be proved as a consequence of Theorem 112.21 
below. □ 

4.3. Estimates for 0. We now turn to the heart of the matter, which is the covariant wave 
equation ( fl8]) for 0( m+1 ). This cannot be viewed as an equation of type □</> = perturbative, 
and therein lies the difficulty To address this issue we identify the nonperturbative part, 
and add it to the main operator □ to obtain a paradifferential type magnetic d'Alembertian. 
Precisely, we define the leading order paradifferential approximation to the covariant O a 
equation using only the free part v4- rce of A[ m \ which is independent of n: 

(43) U\ = D-2iJ2 P <k -cA^P k d 3 , 

k 

Fortunately, this operator does not depend on m. Then the equation for 0( m+1 ) takes the 
form 



D P0(-+i) 



(m+l) > 



(44) 

where M. (A^ m \ 0( m+1 )) is given by 

M(4( m ^KD) = 2i(A {m) ' a d a <f) {m+1) -J^^k-cA^djPk 

k 

(45) _ | i9ii W^Ki) + i (m),a i H ?i ( m+ i)j 

:= M^A^, (m+1) ) + M 2 (A {m \ (m+1) ) 

We further write 

M 1 = A' + A'o 

where M contains the terms of the form A^dj(p and Ao contains the terms of the form A°do<j). 
We remark that M exhibits a null structure. Indeed, using the divergence free character of 
A; we write: 



A, 



VjAi , 



n 



so that: 

M{A x ,<t>) = A l d l <p = ^.\;(V,A %,</>) . 

i<j 

From here on we take the last expression as the definition of N '. This will allow us to retain 
the null form while discarding the divergence free condition on finer decompositions of A x . 

Now the idea is to first produce a parametrix for n A , and then to estimate the right hand 
side of (|44p perturbatively. The linear estimate for the magnetic wave operator O a is one of 
the key points of the paper, and has the following form: 

Theorem 4.3 (Linear estimates for <ft). Let D A be the paradifferential gauge-covariant wave 
operator defined on line (l4"3"|) . and suppose that nA^ ree = with || A-^ ree [0] ||#i xi 2 ^ e. If e 
is sufficiently small then we have: 

(46) H\\ s i < H0[O]b lxL2 + ||n p A < 



i 



Section [B] is devoted to the proof of this result. 

In order to solve the equation 044|) it remains to estimate the the right hand side of 

(47) \\ M ( Aim) ^ im+1) )\\ N nLHH-h % e *ll^ (m+1) H^ 

which is applied with = (fi( m \ In order to estimate the difference 0( m+1 ) — 0( m ) we need in 
addition to show that 

(48) ||-M(^ (m) ,0) -M{A^\<j>)\\ NnL2{il _ h) < (Cer-'Ms 

which is then applied to = (j)( m \ To prove (147|) and (148!) we peel off some easier cases 
before we arrive at the heart of the matter. 

Step l:(The M 2 term.) For this it suffices to have the following estimates: 

(49) II^^W^nL^) ~ ll^oll^^^^ll^lU! , 

(50) II AM) 2 0lh i2(i j-^ )ni i (L2) % II I' A o } 11^^1^00(^1)11^1151 • 

1=1,2 

which are proved using only global Strichartz type bounds and Sobolev embeddings. 

Step 2:(High modulation bounds for Ai 1 ) To establish the L 2 (H~z) bound for Ai l we use 

f )27|) and the following two estimates: 

(51) J2w p <^A f rv x p k( p\\ h < 114^11^11011^, 



k 



(52) II^II^A-i) S Mo 11^11011*- 

It remains to prove the N bounds for A4 1 . At first we separately consider M and Mq. 
Step 3:(Peeling off the good parts of M) The first step will be to show that we can restrict 
our attention to low x high interactions in N ' . For this we define the component J\f lowhl by 

k 

To estimate the difference M — J\f lowhl we will prove the null form estimate 

(53) \\U(A x ,4>)-U l ™ h \A x ,(i>)\\ N < HA, Ml 011* • 
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We remark that once this is done, the free part A x ee has been taken care of and will no 
longer appear. 

For the low-high interactions in TV we still have the null structure to use, but the balance 
of frequencies is unfavourable. For the most part we can still bound J\f lowht via a null form 
bilinear S 1 x S* 1 estimate; the exception to this is the case of high x low — > low modulation 
interactions. We group these in an expression denoted by H* J\f lowht } which is given by 

H*M lowhi {A j ^) = £ Hl,Af(Aj,P k ^) , 

k'<k-C 

For the difference we will prove the bilinear null form estimate 

(54) \\N l ™ h \A x ,<t>)-U*U l ™ h \A x ,<t>)\\ N < IIA^MI^Ik , 

in the last section of the paper. 

The high x low — > low modulation interactions contained in %*j\f lowhl can no longer be 
estimated using the S 1 norm of A x , instead we need the stronger Z hyp norm. This leads us 
to introduce the expression 

H*M lowhi (HA x , 0) 

In view of the estimates fT2B|) , (1221) and (|2ip , the N bound for the expression 

U*M lowhi {A x m \<p) - WN^iHA^-V^), 
as well as the corresponding differences, is a consequence of the following 

(55) || U*M lowhi (A x ,<f>) \\ N < || A x \\ zhw \\ . 

After peeling all the good contributions, the only part of Af which has not been estimated 
so far is 

U*N lowhl (UA x m \<p {m+1) ) 

Step 4:(Peeling off the good parts ofAf ) The arguments here follow the same lines as before. 
However, since A solves an elliptic equation, a larger portion of all cases can be treated in 
a direct fashion via Strichartz estimates and Sobolev embeddings, which leaves less in terms 
of bilinear estimates to be proved later. Recall that Af (A ,(f)) = A d t (fi. Modifying slightly 
the setup before we set 

Af^ h \A ,^= Ys N ^ P <*-cQ<k-cAo,P k <l>) , 
k 

Then the difference is estimated via 

(56) Wo(AoA) -K° whi (Ao,<j>)\\L^ < \\A \\y4<P\\s 

The extra step we can take in the case of Ao is to replace A^ by A^' hh . By (|39l) and 
the second part of (|34l) the difference is estimated using 

(57) ||A^ o "(Ao,0) -^"(^,0)11^ < - A h Q %i L i L ~\\<f>\\ s 
The last two steps are similar to the case of A x . First we introduce 

U*N l r h \AfA) = £ T^A/W, P fe 0) , 

k,k': 
k'<k-C 
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In view of (|39|) and the first part of (134}) . the corresponding differences are estimated using 

(58) \\K wh \Af^)-WK wh \A^\<P)\U < iKW^^UWs^ , 

Finally, replace this by H*M l Q Owhi {UA ,(j)). By @2D and the third part of (JM]) for the 
difference it suffices to have the estimate 

(59) || H*K owM (A h \ <P) - n*K owM (HA , <j>)\\ N < || A^ hh - UA^ \\ zM \\ 0« . 

After peeling all the good contributions, the only part of Af which has not been estimated 
so far is 

H*^ owhl (HA { ^\(f) (rn+1) ) 

Step ^.{Reduction of the remaining terms to quadrilinear null form bounds) So far we have 
one portion of M and one portion of Ao left to estimate. We collect the two together in the 
expression 

R (n+i) = _^*A^°"* i (^A( m) ,0 (m+1) ) +WN lowhi (HA ( ™\<p( m+1) ) 

Here we recall that HA^ and UA^ are bilinear expressions in '. 

The key idea is that there is a cancellation that occurs between the two terms above, 
which is why they need to be treated together rather than separately. The trilinear estimate 
that needs to be proved in this case for the trilinear expression 

R{<f)^,<f)^,(f)^) = -H*N , Q° wh \'HA x ((f) ( - 1 \(f)^), </> (3) ) + H*Af lowhi CHA ((j)^\ </> (2) ), <fi 3 ) 

is 

(60) \\R(^\^,^)\\n < H^ll^ll^ll^ll^ll^ , 
Expanding everything into , and performing some algebraic manipulations^ we have: 

R{^\^\^) = d-Qa-Qs, 

where 

(61) q 1 = n^a^n^d^) ■ d a <p {3) ) , 

(62) Q 2 = ^(A^D-^AW^d^) -8$ (3) ) , 

(63) Q 3 = W,*(A~ 1 n~ l d a d i 'H((f)^di(J) ( - 2 ^) ■ <9 a (3) ) . 

For each of these we prove (160]) separately, which concludes our estimates for Jiii uad , 

4.4. Proof of the LP product estimates. Before continuing with the main thrust of the 
paper, we pause here to dispense with the easiest cases of the multilinear estimates listed 
above. These are ([27]), O, flU}, ([38D, ([39]), flU}, (JUJ), (JEQ), ([52]), ([57]), flU, and (J5QJ. For 
the most part these can be broken down into the following estimates, which are immediate 
consequences of Holder's and Bernstein's inequalitites: 

Lemma 4.4 (Core generic product estimates). One has the following dyadic bounds: 

(64) || P k (A kl( p k2 ) || Ll(i2) < 2^—fr»2- 5 l fc -^l|| A kl \\ L2{ . h) \\ <p k2 || L2(W) , 

(65) ii PMiys) w L2{H - h) < 2 5{fc - max ^h-^\\\ 0« ii^^jii $ n L2(U/6>i) , 



4 See the Appendix for a calculation done without the clutter of abstract index notation. 
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Proof of estimate ( )27|) . This follows at once from summing over (1651) and the inclusion: 

(66) S 1 C L 2 (W^) . 

□ 

Proof of estimate ( j28j) . This follows from the inclusion S* 1 • S* 1 C L°°(L 2 ), a simple conse- 
quence of energy estimates and H 1 C L 4 , and then summing over (1651) . □ 

Proof of estimate ( )29|) . This follows from summing over (]64|) after using ( 166]) and the bilinear 
embedding: 

(67) B-BCL 2 {H^), B = L 2 {W 6 ^) nL 00 ^ 1 ) . 

This last estimate is a straightforward application of trichotomy, putting the high frequency 
term in the energy norm. The standard details are left to the reader. □ 

Proof of estimate ( |38|) . This is immediate from Sobolev embeddings involving L 4 ,Ll,L 2 . 

□ 

Proof of estimates (1591) . (HOI) . (I5"lj) . and (I51Z1) . These follow from (JOB]) as in estimates f[27|) 
and (1251) above. We use: 



(68) L 2 (#5) c L 2 (iy 6 '6) , 

as a replacement for (16"6"|) when necessary. □ 

Proof of estimate ( HTj) . First notice that the desired bound for the cubic terms follows from 
(l£3|) and estimates similar to those above, and A _1 £ 1 L 1 (L 2 ) C £ 1 L 1 (L°°). On the other 
hand, for the quadratic term one has: 

II Wfc^) 11^-) < 2^2-*l^-^«|| l|^ 2CT ^ e , 4) || 1U 2C ^^) • 

when k = max{^} + 0(1). □ 
Proof of estimate ( H9l) . This is again an immediate consequence of ( 1641) , ( )65l) , and ( 1661) □ 



Proof of estimate ( )50|) . For the L 2 (H~^) we can use (165]) once we place the product of A 
in L°°(L 2 ) via if 1 C L 4 . On the other hand the L l (L 2 ) estimate comes from ( 16^1) and using 
(16"T|) for the produce of A and (1651) for 0. □ 

Proof of estimate (1561) . First decompose: 

Aodt^-N^^Ao,^) = T 1+ T 2 , 

where: 

(69) T 2 = P <k _ c Q>k-cA d t (j )k 

We'll show each of T{ G L 1 (L 2 ). For Ti this follows from ( 164|) and using: 

II P> fe -cAo || ^ < E 2 " fc 'H P ^oll L2(H ! r IW*ll i2( W) S 2 1&ll*- 

fe' 
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The bound for T2 is similar except we use: 

\\ p <k-oQ>k-cA \\ i;t ^ ) < 2- fc ||9 t A || J . 2( ^i ) . 
The bound for Hallux 2 follows from (I4ip . placing the second factor in L°°(L 2 ). □ 

Proof of estimate (I57I) . This follows from (|4T|) by placing again the second factor in L°°(L 2 ). 

□ 

We remark that at this point we have reduced the proof of our main Theorem [1] to the 
demonstration of the linear estimates in Theorem 14.31 in Section and then proving the 
multilinear estimates (j3l|), (j32||, (P|h (]54]h (|55|), (|5g]h (155]), and finally ([60]) for each of 
the expressions fl6TT) - fl63l) . These bounds are dealt with in the last section of the paper. 

5. Higher regularity and continuous dependence 

In the previous section we have constructed solutions which are small in S 1 for all small 
energy data. Here we complete the proof of Theorem [T] by discusing the remaining issues, 
namely higher regularity and continuous dependence. The ideas here are fairly standard, 
and in particular closely resemble the similar proof for wave maps, see [30] . For that reason 
we merely outline the arguments that follow. 

5.1. Higher regularity. The goal here is to establish the bound 

(70) IK° nlin H^ + Mh» + Pollv- < \UMA[0])\\hn x h^ 

for all small energy data MKG solutions. This is done by a standard iterative procedure, and 
we explain here how to obtain bounds on V X A X , V x Aq, V x 0. In fact, this is accomplished 
by using the same estimates as before. Commencing with A x , A , recall that we have 

whence we see that y^4* lonIin >( m + 1 ) satisfies a wave equation whose source term is a multilinear 
expression like the one for D ^ onlin ' (m+1) but with one factor 0( m ) or Ax replaced by V x <f>^ m \ 
respectively V ' xA^. Using the same multilinear estimates as before, one then concludes that 

.nonlin^m+l) j-y A nonlin,(m) 1 1 , ^ % m 

||VxA — V x n i nils' 1 ~ l Ue *J 7 

||(1 -U)V x Ar Vm ' {m+1) - (1 -U)V x Ar nn ' {m) \\ zhvP < (Ce*r 
with implicit constant depending on || Va;0[O, •] ||^-i xL2 + || V x Aj[0] ||^i xL2 . A similar argument 
applies to V X A , V x d t A . It remains to bound V x (j)^ m+1 \ Here we recall (JHJ) 

D P0(m+l) = M{A^ m \(t) {m+1) ), 

which leads to 

□^(V x (m+1) ) = V x M(A( m \^ m+ V) + 2iY J ^ x P<k-cA Sree ' J P k d 3 ct> {m+1) 

k 

From (11311) . we find with k\ < k<i — C 

\\W x P kl A^d ] P k2 (j ) ^\\ N < 2^- k ^\\P kl A^ e ^\ s 4V x P k J m+1) \\s^ 
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from which we easily infer 

\\2iJ2^P<k-cA free ' j P k d^ m+1 % < UAMsxIIV^IIs! 

k 

On the other hand, the term 

VJW(A (m) ,0 (m+1) ) 

is again a multilinear expression in 1}, with at most 

one derivative factor in each monomial. Thus the multilinear estimates from above apply. 

5.2. Frequency envelope bounds. The S 1 bounds for the small data solutions can be sup- 
plemented with corresponding frequency envelope bounds in a standard manner. Precisely, 
suppose that {c k } is a H 1 x L 2 frequency envelope for the data, in the sense that 

(71) ||*M4.[0W[0])II*i x l» < c k , \\(A x [0im)\\HixL* « llPfcllp < e « 1 
Assume also that {c k } is slowly varying, 

\c k /cj\ < 2^- fc l, 5<1 

Then we have a similar bound for the solutions, 

(72) IM onlin |k < 4, \\P k <f>\\ s i + \\PkA \\ Y i < c k 

The proof of these bounds is a straightforward consequence of the estimates in the previous 
section, since we have off-diagonal decay in all of our multilinear estimates. 

5.3. Weak Lipschitz dependence on data. While establishing Lipschitz dependence on 
data in the energy norm would be desirable, this does not seem to hold because of the free 
wave component of the magnetic wave equation. Instead, we have a weaker bound for the 
difference of two solutions (Ax , 4>^) and (A x 2 \(j)^) J of the form 

(73) 

for small positive 5. 

This is equivalent to a similar bound for the linearized equation. All the components of 
the nonlinearity for which we have off-diagonal decay in the multilinear estimates cause no 
difficulty, and impose no sign restriction on 5. The only difficulty arises in the paradifferential 
magnetic wave equation 

□^0 = 

Denoting by (B, ip) the corresponding linearized variables, we are led to consider the equation 

k 

where B is an H 1_s free wave. By the bilinear estimate (1541) . for the right hand side we have 
the bound 

\\J2^,<kdMN-s < \\B[0]\\^- SX H- S \\<j>\\sh 5>0 
k 

which leads to the desired linearized bound 

< UMWm-sxH-s + \\B[0]\\n- txA - t U\\^ 5 > 
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5.4. Approximation by smooth solutions. Consider a small energy Coulomb initial data 
(A x ,ip)[0] and its regularizations (Ai , ^ m ))[0] = P <m (A x ,ift)[Q]. Denote by (A x ,(j)), respec- 
tively (A x m \ip^) the corresponding solutions. Our aim here is to prove the following: 

Lemma 2. Let {c&} be a slowly varying frequency envelope for (A x , ip) [0] in the energy norm 
H 1 x L 2 . Then 

(74) \U^-A x ^-n\%<Y, c l 

k>m 

The proof of the lemma is straightforward, by using f J73|) to bound the fequencies less 
tham m for the difference, and by using ( 172]) for the higher frequencies of each of the terms. 

5.5. Continuous dependence. Given a convergent sequence of small data (A x , ip k )[0] — > 
(A x , ip)[0] in H 1 x L 2 , we consider their frequency envelopes {c k,m } respectively c k in I 2 . Due 
to the above convergence, it follows that 

(75) lim \c k ' m \ 2 = 0, uniformly in k 

mo->oo ' * 
m>mo 

We consider the corresponding regularized data (Ax ,ip k '^)[0] and the corresponding 
solutions (A fc,(m) ,^ fc ' (m) ). By the relation ([75]) and (J7JJ) it follows that 

( A k,(m) ^k,(m)j _^ ^ ^ in uniformly in k 

On the other hand by the well-posedness theory for smooth data we have the convergence 

Combining the two we obtain the desired local in time convergence 

(A k ^ k )^(A^) mSl c . 

6. Renormalization 
In what follows we work with the selfadjoint paradifferential covariant wave operator 

(76) ^ A = D-2zJ2A j <k ^ c d,P k 

k 

Here V x ■ A = are Coulomb gauge potentials which solve the free wave equation OA = 0. 
Our first goal is to prove estimates and construct a parametrix for the frequency localized 
evolution 

(77) = (0(O),<MO)) = Qa) 

with all functions 0, /, g, h localized at frequency 1. 

In general A ^ 0, so one cannot in addition have dA = 0; in other words the derivative 
interaction cannot be removed via a physical space gauge transformation. The situation 
changes drastically if one views the gauge potentials as pseudo-differential operators. This 
stems from the fact that when viewed microlocally all connections have zero curvature, 
because they contain only one component for each fixed direction in phase space. To write 
the gauge interaction in terms of a potential we need to solve: 

A-i = -TjP t + ^x- 
18 



This is impossible, but if we further restrict ourselves to the region where r ±|£|, we are 
left with the exact solution: 

A-C = #±-(=F|£l,0 > Mt,z,Z) = -L^A-l(A-u). 
Here ip± are real valued and 

UJ = -^-y Ll = d t ±u-V x , A^ = A - (oj ■ V,) 2 , 

Also we have used the relation L^_L^_ = □ + A w ±. As defined ^ + is associated to the upper 
cone {r = |£|} and if)- is associated to the lower cone {r = — 1£|}. 

Quantizing this, one can write the reduced covariant wave equation approximately as: 

D P A « □ - 2i{d a if)±){t,x,D)d a , ±r>0. 

which suggests that one should be be able to remove the gauge interaction through the 
pseudodifferential conjugation: 

D P A « e~ i ^ ± (t, x, D) □e -i ^ ± (D, y, s) ±r>0. 

Applying this algorithm directly does not work well for two reasons. First, the symbol we 
obtain is not localized at frequency <C 1. Secondly, the symbol \P defined above is too singular 
due to degeneracy of the operator A~l ; this corresponds exactly to parallel frequencies in A 
and cf>. 

To remedy the latter issue we take advantage of the fact that in the bilinear null form 
estimates there is a small angle gain. This allows us to tightly cut off the small angle 
interactions between A and <fi in the construction of ^, though not uniformly with respect 
to the A frequencies. Precisely, we define the dyadic portions of if) by 

(78) ^,±(^,0 = -L%A~i(U- >Sk A k -cu) 

Then the full ^'s are defined by 

For the renormalization we will use frequency localized versions of e namely the opera- 
tors 

e^(t,x,D), e***(D,y,a) 

Here we use P(x,D) for the left quantization, and P(D,y) for the right quantization. Also 
the subscript < stands for the space-time frequency localization at frequencies <C 1. 

Our main goal now is to show that the renormalizations are compatible with the S and 
iV spaces. Below we use the notation 

D A <0 : = D " U^-cdjPo 

and analogously for D A . 
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Theorem 3. The frequency localized renormalization operators have the following mapping 
properties with Z £ {N , L 2 , Nq}: 



(79) ef^(t,x,D): Z^Z, 

(80) dte^^x.D): Z^Z, 

(81) e^' ± (t,x,D)e% ± (D,y,s)-I: Z eZ, 

(82) e-^(t,x, £>)□ - □ p A<o e^ ± (t,x, D) : 7V* ± -> eiV 0i± . 

(83) e3*±(t,x,D): Sj S , 



The above theorem allows us to construct an approximate solution for (|76|) as follows: 

0<w = \ E e <°^ (*» ^7^1^ (A 2/, 0)(|I% ± fc) 

(84) ± 

+ e-; fe (i,^)-^*(A!/ lS )/ 

where 

= / e ±i( *- s)|D| /(s)rfs 
io 

represents the solution to 

(atTip|)u = /, u(0) = 
Here if we drop the e ±l< ^ ± operators we simply have the expression of the exact solution for 
the constant coefficient wave equation. Thus the idea is that these operators approximately 
conjugate the covariant wave flow to the constant coefficient wave flow. 

Theorem 4. Assume that f,g, h are localized at frequency 1, and also that f is localized at 
modulation < 1. Then (p app is an approximate solution for n^ <0 </> = /, 0(0) = (f,g), in the 
sense that 

(85) Ha P p\\s <\\f\\N + \\g\\L* + \\h\\ L 2 

and 

(86) Ua P p[0] - (g,h)\\ip + \\O p A<o <p app - f\\ No < + \\g\\ L 2 + \\h\\ L 2 

We remark that <p app constructed above has the same localization at frequency 1 and 
modulation < 1. We also remark that the actual parametrix construction only requires the 
estimates fjT9"l) - fl8"31) . but yields a weaker form of (I8"5"l) with So replaced by N£. Then the 
estimate (1831) serves to provide the additional S$ regularity. An easy consequence of this is 
the following 

Theorem 5. For all f £ N and (g, h) £ H 1 x L 2 the solution to the paradifferential covariant 
wave equation (ITT)) is defined globally and satisfies 

(87) UWgx < ll/ll^^l + \\9\\hi + \\h\\» 

Again, using only (179 1) - (1831) suffices but yields a weaker form of (1851) with S replaced by N*. 

For the remainder of the section we use Theorem [3] to prove Theorems HI EJ The proof of 
Theorem [3] is completed in Sections [HJ [HI after some preliminaries in Section [7J 
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Proof of Theorem [^} The estimate (!85|) follows directly by concatenating (179|) and ( )83|) . 
It remains to consider (I86I) . We begin with the initial data. For the position we have 

QPP (O) -9=W (0, x, D) (D, y, 0) ± /i) - g 



± 



± 



e^ t (0,x,D) ] l,e^(Al/ ) 0)- 1 l [ 



(|D|<7±/0 



and we can apply the L? version of ((HI 
For the velocity we have 



d4a V M -h=\Y, ±e <^^ x ' D )^o(D, y,0)(\D\g±h)-h 



+ [fte^KO, x, I?) e^= (D, y, 0)(|Z% ± ft)- 



f e-? ± (0,x,D)^[^ ± ](Az/,0)(| J D|^±/ i ) 



1 



±e 



-itp± , 
<0 



(0 I i I D)e^(Ay,0)/(0) 



The first line is rewritten as 
1 



o E e <o^ (0. ^) e <1f (A 2/, 0) - 1 (±|Z% + h) 



and then we can use flHTj) . For the second and third lines we use (1801) . Finally, for the last 
line we use (18 ip twice, along with the bound 

ll/(0)||i. < ||/lk 

derived by Bernstein's inequality due to the unit modulation localization of /. 
Lastly, we consider the error estimate. We have 



'app 



O p A<o e-^(t,x,D)-e^(t,x,D)n}4> ± 





± 


D\ 




D 





£P5 



A< e <0 



e^{D,y,a)f-f 
(t : x : D)-e^ 1± (t,x : D)n]<j )± 



+ -K*(i, I ,C ± (D l!/ , s )-i]/ 



D 



D 



where 



+ e^(t,x,D)— l [d t ^](D : y : s))]f] 



^ = W\ [ e± ' t|D ' e <o ± ( D » °) Q D \9±h) + K± ^ (D, y, a)] 



-±±iip± , 
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For the first line we bound <j>± in N± via ( I79|) and □ 1 estimates, and then use the conjugation 
bound ( )82|) . For the second and third we use (!HT!) in N. Finally for the fourth line we use 
dSDH in N. □ 



Proof of Theorem Consider first the frequency localized problems (with k referring to 
spatial frequency) 

^A^k = ft,, 0fc(O) = (g k , h), k e Z 

We solve each of these approximately and re-assemble the solutions 4>ap P ,k to the full ap- 
proximate solution <j) app = J2k $app,k- We cannot immediately apply the preceding theorem, 
since we make no further assumption on the space-time Fourier support of the source /. To 
remedy this, we split 

f _ fhyp | fell 

where f k vp is supported in the region ||r| — |£|| < 2 fe .Then we solve the two problems 

□P J.1 _ fhyp n P jJ2 _ fell 

approximately, the first by using the previous theorem, the second by neglecting the magnetic 
term 2iA 3 <k dj(j)\. We then solve the second equation by 'division by the symbol', i. e. we set 
(from now on k = by scaling invariance and we omit the subscripts) 

where the operator CT 1 is defined by multiplication with ^j^p on the Fourier side. Then 
the bound 

is immediate, and we reduce to solving the problem 

a A <0 <P 1 = f hyp i <P\0) = (g,h)-n^r ll (0) 

which we do by invoking Theorem HI 

It remains to control the additional error generated by our approximate solution 2 , which 

is 



2iA J <0 d j( j) 2 



This is estimated by 



\\^ <0 d 3 ^ LlLlnL2H _, < 114011^11^11^.^x3 « link 

It follows that our approximate solution (j) app = J2 k 4>app,k satisfies the conditions 

\P P A<l>a PP - f\\ NnllL , H -i + 11^(0) " (9,h)\\frx» « WfW Nnl l L 2H-i + WMWfrxV 

Also, observe that the preceding Theorem [4] implies 

ll^appllsi £ WfW N ni^H-h + WfatyWfrxL* 

The proof is now completed by simple iterative application of the preceding to the successive 
errors. □ 
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7. Decomposable Spaces and Some Symbol Bounds 

7.1. Review of the Basic Decomposable Calculus. First we discuss the notion of de- 
composable function spaces and estimates. Recall that a zero homogeneous symbol c(t, x; £) 
is said to be in "decomposable L 9 (L r )" if c = Ee cW , 6 £ 2_N , and: 



c v 1 i a a < oo 



lOn 



(88) £ II 

e 

where: 

< 89 > » C< X( W! ,) = ll(EE-pii^(^ E )' c M iii : 

fc=0 </> 

Here &g (£) denotes a cutoff on a solid angular sector |£|£| _1 — 4>\ ^ # for a fixed G S™" 1 , and 
the sum is taken over a uniformly finitely overlapping collection. We define || b\\DLi{L r ) a s 
the infimum over all sums (I8"8"j) . In [T5] it is shown that the following Holder type inequality 
holds: 

m m 

(90) || JI^ H^(^) S II II 6 > Wwi) , Or 1 ,*- 1 ) = ^GzrVr 1 ) . 

j=l i=l i 

In the sequel we only need a special case of decompositions provided in terms of these norms: 

Lemma 7.1 (Decomposability Lemma). Let A(t,x;D) be any pseudodifferential operator 
with symbol a(t,x;£). Suppose A satisfies the fixed time bound: 

(91) sup || A(t,x;D) \\ L 2^ L 2 < 1 . 

t 

Then for any symbol c(t,x;£) G DL q (L r ) one has the space-time bounds: 

|| (ac)(t,X]D) | \ L U L 2 _KL<J2(Z/-2) <|| C ||z)i?(Lr) , 

(92) 111 111 

— + - = —, o + - = — ' 1 < 9i,92, 9,r,r 2 < oo 
q\ q q-i 2 r r 2 

Proof. Due to the Z 1 summation over 6* in flHHl) it suffices to consider the case c = eg for a 
fixed 6. We further decompose 

C*(t, x,0 = Yl ^ 0, 4(t, x, := &te)c (0) (*> 

By ( 18 9 p each Cg is supported in an angle 8 sector, and it is smooth on the scale of its support. 
Thus by a Fourier series decomposition we can separate variables and represent 

4(t,x : o^J2 d i' j ^ x >i' j (0: 

j>0 

where 

10n 

iK'Mik sr"£u ^ ^) fc cW ii^' i e M ^ 1 

fc=0 
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Due to the rapid decay with respect to j it suffices to consider the contribution to (? e coming 
from a single j, say j = 1. Then eg has the form 



where 

(93) KIIl^ls < 1, \ef\ < 1 
Then we can represent 

(ac){t, x, D)u = d$(t, x) ■ A(t, x, D)e* e {D)u 

The second factor above inherits the L 2 norm from u due to (T9Tj) and the square summation 
in u due to the sector decomposition. Thus 

\\A(t,x,D)ei(D)u\\ L ^ L2x <\\u\\ L2 

The estimate for (ac)(t,x, D)u follows by combining the last bound with fl93|) . □ 

7.2. A Decomposable Calculus for Pseudodifferential Products. In the sequel it will 
also be useful for us treat estimates for products of operators in a modular way. Recall that 
if a(x,£) and b(x, £) are symbols, then a r b r — (ab) r ~ i(d x ad^b) r . This formula is not exact, 
but it leads to an estimate: 

Lemma 7.2 (Decomposable product calculus). Let a(x,£) and b(x,£) be smooth symbols. 
Then: 

(94) || a r b r - (ab) r \\ l ^(l 2 )^li(l 2 ) ^ II (V x a) r ||l'-(l 2 )->^i(l 2 )II ||zw2(£°°) 

where q~ l = Y1p7 • Furthermore, if b = &(£) is a smooth compactly supported multiplier, 
then for any two translation invariant spaces X, Y one has: 

(95) || aW - (ab) r \\ X ^ Y < \\ (W x a) r \\ X ^ Y . 

Proof. To prove the first estimate (1941) we write the kernel K(x, y) of the difference a r b r — (ab) r 
as follows: 

K{x, y) = c n f e^V^'M*, £) (b(y, rj) - b{y, 0)dzd^d v 

= c n [ [ e i{x - zH e i{z - yyri a(z, £)Ve&(y, sr)+(l- s)£) ■ (77 - £)dzd£dr]ds , 

JO Jl!L 3n 

= c n i [ [ e^-^e^-^V^z, ■ V i b(y, sr)+(l- s)£)dzd£dr)ds , 
Jo </R 3 " 

(96) =i 77r^[ I T {s ^V x a T (x,D)T_ sE vlb(-,~)dsdZ, 
where we have used Fourier inversion: 

V € 6(y, sr, + (1 - = t^- / e 4 ^ 1 -^ vS(y, S) <E . 
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From the smoothness of b on the unit scale in £ we have that the weight function obeys 
the estimate: 



/ II V^6(-, H) || LP2 (L°°)rfH < || 6 ||diLP2(l°°) 



Then ( 1941) is a direct consequence of Holder's inequality and the the translation invariance 
of LP spaces. 

Note that the proof of (1931) also follows from the identity (1931) directly because in this 
case the weight function V^fr is a constant in y. □ 

7.3. Some Symbols Bounds for Phases. For its use in the sequel, we list out a num- 
ber of decomposable estimates for the phase ?/>(£, x;£) used to define our microlocal gauge 
transformations: 

Lemma 7.3 (Decomposable estimates for ip). Let the phase ij)(t,x;£) be defined as in (1751) . 
and its angular components ^ = rr;£), where cu = |£| -1 £. Then for q > 2 and 

2/q + 3/r < 1 one /ias: 

(97) || (tf>, 2- fc V f >f) IU W S 2-(i^) fe ^-!-?e, 
In particular 

(98) 11(^,2-^,^)11^^00) < 2^ fc e, g>4, 

(99) || |U W < 2<*-*-'<i+*»* e , r > 6 , 

Proof. Notice that the last two estimates follow from the first by summing over dyadic 
2~sk < g < 1. For the first bound we interchange the t integration and the u summation to 
obtain: 



where at the second step we have used the Coulomb gauge to gain another factor of 9. 

Now we conclude using the Strichartz estimates. In four space dimensions the Strichartz 
sharp range is given by - + ^ = |. Moreover, on an angular sector of size 6 Bernstein's 

0/^_ 1_ \ At 1 j_\t> 

inequality gives the embedding Hg(D)P k L ro C 9 r <> r>2 (r o r> U . Thus: 



k \\S k , 

and the second estimate follows. □ 

We wrap this section up by proving some additional symbol type bounds for the phases 
ip. These involve the variation over the physical space variables: 

Lemma 7.4 (Additional symbol bounds for ip). Let ip be as above. Then one has: 

(100) |V< fc (t, x- - ?P <k (s, y,£)\<e log(l + 2 k (\t - s\ + \x - y\)), 

(101) \i[>(t, x- - 4j{s, y;Z)\< elog(l + |t - s| + |x - y|) 

(102) |<9^(t, s; - i[>(s, y; 0)1 < e<(* - s, a; - y)) 1 ^, i < a ^ a~\ 
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Proof. We decompose as before 

For each fixed # and j we have by the definition of ijj and the Coulomb gauge condition 

(t,x,0l<^ 1 2- J sup||n^|| i oc 

Then by energy estimates for /I and Bernstein's inequality we obtain 

(103) \^ 6) (t,x,t)\ < etWAMWmxv, |^(i,z,OI < ||^[0]||^ lxL2 
A similar argument leads to 

(104) \d t ^f (t,x,0\ < 2^l||A i [0]|| iflxi2 , |c^V>;(M,OI £ ^IIA^Olll^i^ 
Differentiating with respect to £ yields 6 1 " 1 factors, 

\d^f\t,x,0\ < O^WMO}^^, \d x , t d^?\t,x,{)\ < 2W^\\A,[0]\\ mxL2 . 

For the bound (11001) we use both fllQ3j) and (17.31) to write for j < k 

\ip<k(t, x; - ^< fc (s, y- i)\ < 2 j {\t -s\ + \x- y\) + \k - j\ 

and then optimize the choice of j. 
The proof of (I102p is similar. 

□ 

8. L? ESTIMATES FOR THE GAUGE TRANSFORMATIONS 

In this section we prove three core L 2 based estimates for the gauge transformations e ±l ^ ± . 
These will serve as building blocks in later sections. 

8.1. Oscillatory integral estimates. In order to prove various estimates involving the 
operators e<g (£, x, D) and e<Q(D,y,s) we need to obtain pointwise kernel bounds for 
operators of the form 

T a = e -^±( t ,x,D)a(D)e ±t{t - s)lDl e^ ± (D,y,s) 

where a is localized at frequency 1. The kernel of the operator T a is given by the oscillatory 
integral 

K a (t,x;s,y) = j e^(t,x^)a(^^e^^^,y,8)d^ 
Our main estimates for such kernels are as follows: 

Proposition 6. a) Assume that a is a smooth bump on the unit scale. Then the kernel K a 
satisfies 

(105) \K a (t,x; s,y)\ < (t - S )-% - s\ - \x - y\)~ N 

b) Let a = ac be a bump function on a rectangular region C of size 2 k x (2 k+l ) 3 with 
k < I < 0. Then 

(106) \K a (t,x; s,y)\ < 2 ^ l (2^ k+l \t - s))^(2 k {\t - s\ - \x - y\)y N 
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If in addition x — y and C have a 2 k+l angular separation then 

(107) \K a (t, x; s, y)\ < 2 ^ l (2 2 ^\t - s\)- N {2 k (\t - s\ - \x - y\))- N 

Proof, a) Away from a conic neighbourhood of the cone {\t — s|=±|x — y\} the phase 

\f = ±(t - s) |e| + £(x -y)- (Mt, x, o - y, 0) 

is nondegenerate due to (j!02p with \a\ = 1. Hence repeated integration by parts yields 

1^, s,j/)| < ((*, x) - (s, y))~ N , N ~ a" 1 

Near the cone we need to be more careful. Denoting T= |t— + y| and i? = \t— s\ — \x— y\ 
in suitable (polar) coordinates this takes the form 

K a (t,x jSl y)= f e -«(^(*.^0-^(*.i/.O) e iB& e fl'€' 3 a ^)^ 



In (the former radial variable) this is a straight Fourier transform. Given the bound fU02p . 
we can use stationary phase in While the £ derivatives of the ^± part of the phase are not 
bounded, they only bring factors of T a , which is small enough not to affect the stationary 
phase ( this works up to a — |). We obtain 

\K a (t,x,s,y)\<T-l(l + R)- N 

b) Away from the cone the estimate follows easily as above since the phase is nondegen- 
erate. Near the cone we use again polar coordinates to express our oscillatory integral as 
above, 

K c (t,x,s,y) = [ e -*(^(*.*.0-^±(-,»,0) e iBSi e «'€' 2 ac( ^ t ^ 



where ac is a bump function in a rectangle on the 2 k scale in the radial variable £i and on 
the 2 k+l scale in the angular variable Then we can separate variables in We note 

that this rectangle need not be centered at £' = 0, though this is the worst case. In £i this 
is again a Fourier transform, so we get the factor 

2 k (2 k R)- N 

In £' we can use stationary phase to get the factor 

23(fe+0/2 2 (fe+0 T \-| 

The bound (11061) follows by multiplying these two factors. 

Finally, the estimate fjlOTj) corresponds to the case when ac is supported in > 2 Z in the 
above representation. If T < 2~ 2 ( fc+/ ' ) then there are no oscillations in £' on the 2 fc+/ scale, 
and we just use the brute force estimate. For T > 2~ 2 ( fc+/ ) the phase is nonstationary in 
and we obtain the factor 

□ 

While the above proposition contains all the oscillatory integral estimates which are 
needed, it does not apply directly to the frequency localized operators (t, x, D) and 
e^ ± (-D, y, s). For that we need to produce similar estimates for the kernels K a <0 of the 
operators 

T< = e2 i ' ± (t,x,D)a(D) e ^ t - s ^ e ^(D,y,s) 
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The transition to such operators is made in the next 

Proposition 7. a) Assume that a is a smooth bump on the unit scale. Then the kernel K^ 
satisfies 

(108) \K (t,x; 8 ,y)\ < (t - S )~% - s\ - \x - y\)~ N 
In addition, the following fixed time bound holds: 

(109) \K a <0 (t, x) t, y) -a{x-y)\< e\ log e| 

b) Let a = ac be a bump function on a rectangular region C of size 2 k x (2 k+l ) 3 with 
k < I < 0. Then 

(110) \K%{t,x; s,y)\ < 2 ^ l (2^ k+l \t - s))^{2 k {\t -s\-\x- y\))~ N 

c) Let a = ac be a bump function on a rectangular region C of size 1 x (2') 3 with I < 0. 
Let wG§ 3 be at angle I from C . Then we have the characteristic kernel bound 

\K%{t, x; s, y)\ < 2 3 <<2 2 <|t - s|)-"<2y - yf\)~ N 
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Proof, a) We represent the symbol as 



(x-y)-u 



112) e<? ± = J m{z)e ±iT ^ dz 



where m(z) is an integrable bump function on the unit scale and T z denotes translation in the 
direction z, with z representing space-time coordinates. Since the wave equation is invariant 
to translations, the symbol e ±tTz ^ ± is of the same type as e * ,± . Using this representation 
for both ijj± exponentials, the kernel _fQ can be expressed in the form 

K%(t, x,s,y)= [ [ e~ iT ^± (t, x, £) a (£)e l(±m >(t-^-y) e ^± ( S) y? ^ ^ m (z)m{w) dzdw 
T Z T W [ e- iqp± (t,x^)a(^)e l{±mHz - w) e l ' p± (s,y^) d£ m(z)m(w) dzdw 



Denoting a(z,w)(£) = a(^)e^ ± ^' , ^- ) '^ w \ we can express the above kernel in terms of the 
kernels K a in the previous proposition, namely 

(113) K^ (t,x,s,y) = J T z T w K a{z ' w \t,x,s,y) m{z)m(w) dzdw 

To prove the bound (11081) we use (11051) . together with the additional observation that the 
implicit constant in (1105)) depends on finitely many seminorms of a (at most 8, to be precise) 
which we denote by 1 1 1 a \ \ \ . Then 

\\\a(z,w)\\\<(l + \z\ + \w\) N 

However, this growth is compensated by the rapid decay of m, therefore the bound (1105)) for 
K a transfers directly to if< in (11081) . 

To prove (11091) we use the same representation as above to write 

K%{t, x, t, y)-a(x-y) = J J [ e -<(^±(*,*,0-T»iMt.v,fl) _ i] a ^)e^ x ~^ d£ m{z)m{w) dzdw 
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By ffTOTj) we have 

\T z ip±(t,x,£,) - T w ip±(t,y,£,))\ < elog(l + \z\ + \w\ + \x - y\) 

which yields 

\K^. (t, x, t, y) — a(x — y)\ < e J log(l + \z\ + |iy| + |x — y|)|m(2)||m(u;)|(i2;(iw 

< elog(2 + \x- y\) 

This suffices if log(2 + \x — y\) < | loge|. But for larger |x — y\ we can use (11081) directly. 

b) Using the representation (I113p . the bound (11 101) follows from H106[) exactly by the same 
argument as in case (a). 

c) Using the representation f 1 1 1 3 j) . the same argument also yields the bound (11051) provided 
we have the following estimate for K a : 

\K a (t, x, s, y)\ < 2 3l (2 2l \t - s\)- N (2 l \x' - y'\)- N (l + \(t - s ) - (x - y) ■ co\) 10N 

To see that this is true, we consider three cases: 

(i) If \t - s\ < 2~ 21 then fTl06D applies directly. 

(ii) If \t-s\ > 2~ 21 but \\x-y\-\t- s\\> 2 l \x' -y'\ + 2 2l \t - s\ then (fT06D still suffices. 

(iii) If \t - s\ > 2~ 21 and \(t - s) - (x - y) ■ cu\)\ > 2 l \x' - y'\ + 2 2l \t - s\ then fTl06D also 
applies. 

(iv) Finally, if \t - s\ > 2~ 21 , but \\x - y\ - \t - s\\ < 2 l \x' - y'\ + 2 2l \t - s\ and \(t - s) - 
(x — y) • uo\)\ 2 l \x' — y'\ + 2 2l \t — s\ then we must have Z(x — y,u) 2 l , which implies 
that Z(x - y, C) w 2 l . Then fTTOTj) applies. 

□ 

8.2. Fixed-time L 2 estimates. The following is an important application of the previous 
theorem which is at the heart of our parametrix construction: it proves the L 2 -part of (17U|) . 
PD as well as that of (|8TD. 

Proposition 8. The following fixed time L 2 estimates hold for functions localized at fre- 
quency 1: 

(114) e**±(t,x,D): L 2 -+ L 2 , 

(115) e£±(t,x,D)e$t{D,y,8)-I: L 2 ^e^\ogeL 2 

(116) d x ,ef^(t,x,D): L 2 -+ L 2 
Proof, a) By the estimate (I105P with s = t, the TT* type operator 

e ± ^(t,x,D)P 2 e^(D,y,t) 

has an integrable kernel, so it is L 2 bounded. Therefore e ^(t, x, D)P$ and its adjoint are 
L 2 bounded. To accomodate symbol localizations we observe that 

±itp± I t \ ±iT z ip± 



e <k ± = J m k (z)e <0 zW± dz 

where m(z) is an integrable bump function on the 2~ k scale and T z denotes translation 
in the direction z, with z representing space-time coordinates. Since the wave equation is 
invariant to translations, the symbol e ±lTz ^ ± is of the same type as e ± *^ ± and its left and 
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right quantizations are also L 2 bounded. Thus the bound (11141) follows by integration with 
respect to z. 

b) For the estimate (11 1 5 [) we note that the kernel of (t, x, D)a(D)e<Q (D, y, t) — a(D) 
is given by K< Q (t, x, t, y) — a(x — y). Combining fl 105[) and (I109P we get 

\K<o(t,x,t,y) - a(x - y)\ < min{e| loge|, \x - y\~ N } 

JV — 4 

The integral of the expression on the right is about e~^|loge|, therefore the conclusion 
follows. 

c) By translation invariance we discard the < symbol localization, and show that 
dx,te ±l1,J± {t,x, D)P is L 2 bounded. We have 

By flHH]) we have d x ^± G DL 00 (L 00 ) therefore we can dispose of it and use the L 2 bounded- 
ness of e ±i ^ ± (t,x,D)P . 

□ 

8.3. Modulation localized estimates. The next order of business is to show that the fixed 
time L 2 bounds for e^g* drastically improve to space-times L 2 (L 2 ) bounds if one selects a 
fixed frequency in the symbol. Precisely, 

Proposition 9. For I ^ k' ± 0(1) one has the fixed frequency estimate: 
(H7) || Qie$*Q<oPo II o,i < 2^ fc V 

TV*— >X 1 1 

In particular summing over all (I, k') with I ^ k and k — 0(1) ^ k' for a fixed k ^ yields: 
(H8) || Q<k(e±i* - ef* c )Q <0 P || , 4 < e . 

iV — ^ -A ^ 

A key step in the proof of the proposition is the following result, which we state separately 
since it is of independent interest: 

Lemma 10. Assume that 1 ^ q ^ p ^ oo. Then for k + C < I < we have : 

(119) || (ef^< fc )(t, x- D) \\ LP{L ^ Lq{L2) < e2^ k 2 5 ^ , 

This holds for both left and right quantizations. 
Proof. For the symbol we iteratively write: 

Sie *i><H = ± i2 - l S l (d^ <k -e ±i ^) , 

5 

= ... = (±z) 5 2- 5L l[[Sl j) d^ <k }-e ± ^ , 
i=i 

where the product denotes a nested (repeated) application of multiplication by Sid t ip<,k, for 
a series of frequency cutoffs S^^Si = ~ Si with expanding widths. Disposing of 
these translation invariant cutoffs we see that (11191) follows directly from (198j) . 

□ 
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Proof of Proposition 0. We proceed in a series of steps. 

Step l:(High modulation input) First we estimate the contribution of Qk&w Q>k-cPa to 

o i 

line ( I117p . Using the X^ 2 bounds for the input, it suffices to prove the estimate: 

II Qk^Po Wlhl^l^) < 2*<*-*V 
By Sobolev estimates in |r| ± |£|, this reduces to the bound: 

II „±i* o|| <- n—hk' 

\\ e k> P o\\ L 2 {L2) ^ {L2) S 2 o e. 
Using continuous Littlewood-Paley resolutions to decompose the group element we have: 

Jk">k'-C 

(120) = £ + . 

We'll treat these two terms separately. In the case of £ we use estimate (|98|) with p = 5. To 
estimate the contribution of 7?. we use estimate (I119p . 

Step 2:(Main decomposition for low modulation input) Now we estimate the expression 
QkG k ^Q<k-cPoU- First expand the untruncated group elements as follows: 

e ±i * = e ±m <"- c ±i [ ^ t e ±m < k - c dl- [[ V l %,e ±i9 <''- C dldl' 

Jl>k-C J J l,V>k-C 

Tiffl ^i^i^i"e ±m <i"'dldl'dl" , 

JJJ l,l',l">k-C 

= Z+C+Q+C . 
We will estimate the effect of each of these terms separately. 

Step Z\(Estimating the zero order term Z) After localization via S k > the desired estimate 
follows directly from ( I119p . 

Step 4:(Estimating the linear term C) First split the integral as follows: 

(121) C=f ^ l S <k , c e^<^dl+ [ ^> l S >k , c e ± ^-odl . 

Jl>k-C Jl>k-C 

Step 4a:(Estimating the principal linear term in C) For the first term on RHS of line (11211) 
it suffices to show the general estimate: 

(122) || Q k ■ (ipib <k _i c )(t, x; D) ■ Q <k -cPo \\l°°(,l*)-h?(&) 

< e2- 1 > k 2^ k - l hu V \\B <k i c (t)\\ L2 _> L2 , 

t 2 

for / ^ k, and for symbols b(x,C,) <k _i c with either the left or right quantization. In this 

case the modulation of the output determines the angle between the spatial frequencies of 

ipi(x,£) and the spatial frequency of the input, which is 9 ~ 2 2 ( fc— Since this is also the 

angle with £, we may restrict the symbol of to ^2 &>2 ^(k-i) ipi^ for which the estimate ( I122p 

follows immediately from ( 19 2 p and summing over ( 1971) . 

31 



Step 4h:(Estimating the frequency truncation error in C) For the second term on RHS of 
line (I12ip we use (|98|) for ipi with p = 6 combined with fll 19j) with (jp2, q) = (6, 3). If I < k' — C 
the Sk> localization lands on the exponential and we gain 2~2 fc 2e 
we can disregard and directly get 2~2 k 2e( k ~ l \ 

Step 5:(Estimating the quadratic term Q) We follow a similar procedure to Step 4 above. 
First split S <k _iQe ^ <k -° + S >k _i c e ^ <k -° . For the second term one can proceed as in 
Step 4b above using ( 198|) . (11191) . and ( 190|) . Therefore we only need to consider the effect of 
the first term, for which we'll show the trilinear bound: 

(123) || Q k ■ faM>i'&< fc _ic)(*> x > jD ) ' Q<k-cPo IU~(£2)->L2 (i 2) 

< e >2-^iM2iM Bup ||B <Mc (0||x»-, £ . ) 

t 2 

for V ^ I ^ k. By the localization of the output, the symbol ipi^v ls localized to the sum: 

6) > 2 ^(fc-i) 6) /> 2 ^(fc-i') 

For the term Ti put the first factor in DL 3 (L°°) and the second in DL 6 (L°°). This gives us 
dyadic terms in LHS (ll23p (Ti ) ~ 2~^ k 2^ k ~ l ">2^ { - k ~ l '\ For the term T 2 do the opposite, which 
yields a similar bound. 

Step 6:(Estimating the cubic term C) In this case we can gain 2^ k ~ k " 1 directly through the 
use of ( 198]) and three DL 6 (L°°) unless I, I', I" < k' — C. In the latter case Sk> localizes the 
exponential and we can use a ( 11191) instead. Further details are left to the reader. 

□ 



9. Proof of the N -> N bounds in (J79D and (180]) 

Because of use in the sequel, we'll prove a somewhat more general and symmetric version: 

Proposition 9.1 (Symmetric N Estimates). One has the following operator bounds for 
either right or left quantizations: 

(124) ||e±f \\ N ^ N < 1 . 
as well as 

(125) \\d t e±f\\ N ^ N < 1 . 
In particular, by duality one also has: 

(126) || e±f \\ N *^ N * < 1 , 
for either right or left quantizations. 

Proof. Applying e<Q* to F E N we need to consider the following cases: 

Case 1:(F is an L X (L 2 ) atom) This follows at once from the fixed time L? — > L 2 mapping 

property of e^*- 
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0,-i 

Case 2:{F is an X^ 2 atom) Fixing a modulation it suffices to consider Q k PoF, and show 
that: 

\\ef*Q k P F\\ N < 2-i fc ||F|| i2(i2) , 
which by duality reduces to showing: 

II Qu^qPoG I) l 2 (L 2 ) < 2 _ 2 fe || G , 

for a slightly larger cutoff P . At this point the estimate follows the usual game of splitting 
e <o* = e <fc-c + ( e <o" 1 ' — e <fc-c)- For ^ ne ^ rs ^ term there is no modulation interference and 
the estimate is direct. For the second use (11181) . □ 

10. Proof of the Renormalization Error Estimates 

10.1. The N and N* Space Recovery Estimate (EI]). We prove (JED for N*. The L°°L 2 
part of the estimate ( IHTj) for N* is a direct consequence of the similar fixed time L 2 bound 

0,1 

(I115p . It remains to consider the Xoo 2 estimate. Denote 

R k = e^ ± (x,D)e i ^(D,y) 

Fixing the modulation we write: 

Qk(Ro — I) — Qk(Ro — i)Q>k~c + QkRoQ<k-c 

For the first term we bound Q >k _c i n L 2 and use again (j!15p . For the second we can freely 
subtract the low frequencies in the symbol 

QkRoQ<k-C — Qk(Ro — R<k-c)Q <k-C 



Q k (e^(x,D) - e^_ c (x,D))e% ± (D,y)Q <k ^ c 



+ Q k e<7-c& D)(e7o(D, y) - e^_c(D, y))Q<k-c 

o,i 

We estimate the two lines separately as operators from N* to X^ 2 . In the first line we can 
use ( I126p to discard e l ^(y, D)Q <k _c, and then (11181) for the remaining factor. In the second 
line we first write 

Qke~ k %(x, D) = Q k e~' k %(x, D)Q k 
and use (11141) to discard the two left factors. Then we are reduced again to (I118p . 

10.2. Proof of the conjugation estimate (1821) . Recall the following product formula for 
left quantizations: 

a(x,D)b(x,D) = c(x,D), c(x,0 ~ e'^^ (a(x, 0Kv,v)) k= y • 

We are interested in the case where the first factor is the d'Alembertian. Then the above 
formula contains only finitely many terms and is exact and we'll denote it by Using this 
to compute both terms on line (1821) we have: 

i-d 2 t - |e| a )#«e^ =e^{-d 2 t - \i\ 2 ) + 2[id t e^d t - %e^] - Ue$ 

=e<? ± (-^ 2 -|e| 2 ) + 2^e-? ± (^±|e|) 

]<o 



2i[{±d t ^ ± \i\-d x ^)e- i 
[W±\ 2 -\d x ^ ± \ 2 )e-^U 
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On the other hand we can write 

iA^e^ = i[A j e-^U^j + AWje^ + i[A 6 , S^e'^j 
Hence combining all terms we can write down the exact symbol for the conjugation operator 
Diff := (-d? - K| 2 + 2^-)#ze-? ± - e$*(-% - |£| 2 ) 
= - 2z[(±^±|£| - d x ^ + AiQe'^U 
+ 2d t e-^{idt±\t\) 

-2iA j [d^ ± e-^U 

+ 2t[A J ,S <0 ]e~ i ^ 
:= Diff! + Diff 2 + Diff 3 + Diff 4 + Diff 5 

It remains to estimate all five terms as pseudo differential operators from S± C N± to N±. 

The estimate for Diffi: We recall that ip± was chosen exactly so that the high angle 
interaction part of Diffi cancels. Using the definition of ip± we rewrite this term as 

Diff 1 = [(J-n CT )^-e)e- #± ]<o 

We can replace this by 

Diff^p-n^A-Oe-^]^ 
where the inner < truncation is on a slightly larger ball. By translation invariance we can 
discard the outer < truncation. Also here we do not need the S" structure, and it suffices 
to work with N*. Then we can also use the bound ( 1126 h to discard e<g . Hence we are left 
with having to prove a bilinear estimate, namely 

(127) (I - W)A j dj : N* -»■ N 

After a frequency and angular expansion of Aj, we obtain a sum 

fc<0 l<ak ui 

where the multipliers Pf and Pf select 2 l sectors which are 2 l separated. To estimate the 
summands in N* we split the into low and high modulation, 

P?A{pfu, = P?A{PfQ >k+2l _ c d 3 u Q + P?AlPfQ <k+2l _ c d 3 u 

For the first term on the right, we place it into L\L 2 X by using L\U^ for the factor P^A 3 k 
and l\ x for the second factor Q^+21-cdjUo. The Coulomb condition results in a gain 2 l 
(due to the operator dj), and Bernstein's inequality via PkPj^L 2 L x C L\L^ yields a gain 

k-\-l I 

2~ , whence recalling the definition of iV*, we obtain a net gain of 2a, which suffices for 
summation both over I and k due to the condition I < ah. 

■ -i 

For the second term on the right, we place the output into 2 , which is possible since 
its modulation is of size ~ 2 k+21 from elementary geometry. Placing the second input djUo, 
the numerology is the same as for the first term. 

The estimate for DifF 2 : By the definition of S± we have 

id t ± \D\ : S± : N 
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Then we can use directly the bound (11251) . 

The estimates for Diffs and Din^: Here we use decomposability. Precisely, by ( |97|) we 
have dt, x ip± £ DL 2 (L°°). It is also easy to see that A belongs to the same space. Then by 
(190)) we can place the output into the energy space L\L 2 X . 

The estimate for Diff 5 : This term is handled exactly like Diffj; this time we do not have 
the small angular separation condition between the factors, but summation over k becomes 
possible since we essentially gain an extra low-frequency derivative: 

[Aj, S^jdjUo = L(yt,xA\djUQ) 

where L is a translation invariant bilinear operator with integrable kernel. 

11. Proof of the Dispersive Estimates 
In this section we prove the mapping property fl83l) in a series of steps. 

11.1. Proof of the basic Strichartz estimates. First we show: 
(128) e^(t,x,D):S{ ± -> S str . 

The first step is to reduce the problem to the case of homogeneous waves. This is done by 
foliating with respect to t for the L l L 2 part, and with respect to {r — |£| = cr} slices for the 

Xi' 2 part, 



«*(«) = / e ito e ±i, l D lp„/ ± (r)d M ± (r) , / || /*(t) Wwdffr) < 



Then it remains to show that 



e^ ± (t,x,D)e ±it \ D \p : L 2 -> S str 

By a TT* argument this is equivalent to 

e2 i ' ± (t,x,D)e ± ^- s ^P e i ^ ± (D,y,s) : U> V UU 

for all pairs of Strichartz exponents (p,q). In the non-endpoint case p > 2 his follows in 
a standard manner from an L 2 — > L 2 fixed time bound (see (1 1141) ) and and the dispersive 
estimate (I108p . In the endpoint case L 2 L 6 we can use Theorem 4 in [28], which asserts that 
the endpoint case follows from the non-endpoint case plus the dispersive estimate. 

11.2. Proof of the square summed Strichartz estimates. Fixing an angular scale 2 l 
and the corresponding modulation scale 2 21 as well as a rectangle scale 2 k x (2 k +l ) 3 with 
k' < k,l' < and 21 < k' + I' < I we need to show that 



|| J Pc , Q<2ze < o V ' ± (t,x,L')M|| 2 est , < || 



For this we write 

<2<2«e~o V,± (t, x, D)u = Q<2i(e<o V,± - e~<ti-c) u + Q^e^c^^+c^ 

The first term is estimated in X^' 2 using (11181) . For the second one, the symbol of e^i-c 
is strongly localized so that non adiacent rectangles do not interact. Then the square sum- 
mation with respect to cubes is inherited from S\ and it remains to prove that for a single 
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cube C we have 

(129) \\e2tc^x,D)P c u\\ S str < \\u\\ st± 

We can freely discard Pc- Then the last bound follows by translation invariance from (11281) . 

11.3. Proof of the square summed L 2 L°° estimates. The setting and the argument is 
identical to the one above up to the counterpart of (11291) . which now reads 

(130) \\e^ c (t,x,D)Q <2l+c P c u\\ L 2 L o <2 k '\\u\\ s » ± 

We can no longer discard Pc due to the presence of the scale factor 2 k> . Instead we repeat 
the argument in the proof of (j!28p . irs the problem reduces to an estimate for free waves, 

\\e^t c (t,x,D)e ±it ^P c u\\ L 2 Loa < 2 k '\\u \\ L 2 

Then by a TT* argument this is equivalent to 

\\e^(t,z,D)^^WPce%{D,v,8)f\\v^ < 2 2k ' \\f\\^ 

This in turns follows from the dispersive estimate (IllOj) . which shows that kernel of the 
operator above is bounded by 

2 4k ' +31 ' (2^ k ' +l '\t - 

which integrates to 2 2k ' +v < 2 2k ' . 

11.4. Proof of the plane wave norm estimates. Again using the nesting (jSJ), at modula- 
tion < 2 21 it suffices to consider expressions involving Q <2 ie^^ c (t, x, D). By foliating as in 
the proof of (11281) . we further reduce considerations to expressions e~ 2 ^ c (t, x, D)e ±lt ^P f 
for || / 1 1 2,2 ^ 1. Localizing e^_ c .e ± **' £) 'Po/ by Pc, we employ polar coordinates in Fourier 
space to write: 

e^ c (t,x,D)e^P c f = ^J^*!%^ffl,x)<kj, 

where: 

f£(t,x) = Je i{±t+xuj)x p c (\uj)f(\uj)\ 3 d\. 

Then computation of the square sum PW T norms follows immediately from these expressions, 
Holder's inequality, and the Plancherel's theorem. 

11.5. Proof of the null energy estimates. As in the case of the square summed Strichartz 
estimates, the problem reduces to proving that 

e"^ (£,£,£) : S 1 * NE 

We fix a null direction uj and prove the null energy estimate in its associated frame, 

fe-^(t,x,D):S«^L~Ll ± 

For simplicity we take u = (1, 1, 0, 0, 0). The symbol of its associated tangential gradient y 
then consists of r — £i and Depending on whether we are near the + cone or near the — 
cone we can replace 

7--& = (t±|£|)-(6±I£I) 
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The contribution of the first term to the output is easily estimated via the N* norm in L 2 
and then in L^L 2 ± by Bernstein. It remains to bound 

(D 1 ±\D\,D')e^(t,x,D):Sl^L^Ll ± 

We fix the signs to +. Since no time multipliers are involved any more in the above operator, 
the reduction to the case of homogeneous waves still applies, and we are left with proving 
that 

(A ± \D\,D')e-^(t,x,D)^ D \ : L 2 -+ L™L 2 ^ 

The obvious course of action at this point is to use disposability arguments to commute the 
derivatives inside, and then use a TT* argument. This unfortunately borderline fails due to 
the limited decay of the wave kernel along null cones. Instead we take an extra step, and 
perform a dyadic decomposition with respect to the angle between £ and to' = (1,0,0,0). 
There are two cases we need to consider separately. 

For large 0(1) angles we can simply discard the multipliers, and show that 

P?e^(t,x,D)eW:L*^I£Ll. 

For this we commute P^' inside, using disposability to estimate the error in L 2 . Then, by 
translation invariance, it remains to show that 

e^ ± (t,x 7 D)e^Pf : L 2 L 2 (£), £ = {t = 00' ■ x} 

By a TT* argument this is reduced to an estimate of the form 

e^{t,x,D)^Wpi'e^(D,y,s) : L 2 (S) L 2 (£) 

Due to the angular separation, we can use (11071) with I = to conclude that the kernel of 
the above operator decays rapidly off diagonal, and the L 2 boundedness easily follows. 

For small <C 1 angles we can write D\ ± \D\ = a(D)D' and then discard a, and we are 
left with proving that 

D>e-^{t,x,DyW-.L 2 ^L~Ll ± 

The advantage of this is that multipliers which depend only on D' are compatible with the 
L 2 X . Hence via a dyadic decomposition in D', we need to show that 

2 k P k (D')e^ ± (t,x J D)e im : L 2 — >■ /^L 2 (S), k < -C 

We commute the multiplier inside, 

2 k P k (D')e-J ± (t,x,D)e u W = P k (D')[2 k P k (D'),e2^(t,x, D)}e u ^ 

+ P k {D')e^t{t,^D)e a \ D \p k {D') 

By the commutator estimate ( I94p the first term satisfies the same bounds as the operator 

A.(^(V'e<f)(^,i))e lt|D| = -^(Z)0(V^e- #± )<o(^ 1 I))e tt l^ 

Since Vip± G DL 2 L r for some r < 00, by disposability we have 

{V^±e- itp± ) <0 {t,x,D) : L°°L 2 -)• L 2 L^ 

Hence by Bernstein, 

\\hm{Ve^){t,x,D)e^\\ Ii ^ < 2-~r k 
where the point is that there is some gain in k in order to guarantee summability. 
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It remains to show that 

P^Dy-^it^Dy^P^D') : L 2 Z*L 2 (E) 

We harmlessly discard Pk(D'). The square summability with respect to k is now provided 
by the multiplier on the right, so it suffices to fix k. By a TT* argument this reduces to 

e^{t,x,D)e^ D WP*e'g{D,y,8) : L 2 (S) -> L 2 (S) 

For the above operator we use the kernel bound provided by fillip with I = k, namely 

2 5l (2 2l \t-s\)~ N (2 l \x-y\)- N 

This is an integrable kernel, so the L 2 — >■ L 2 bound easily follows. 

12. Statements and proofs of the core multilinear estimates 

It remains to prove (j3Tjh (|32j), (P]l, (155]), (153]), (155]), (f55]> ; (Pjh and finally QSDJ for each 
of the expressions (j6"Tj) - (}53")) . This set of estimates can be roughly grouped into two rough 
categories. The first involves "null-form" type estimates, and the second involves product 
estimates without a microlocal gain from angular separation. The first category of bounds 
boils down to: 

Theorem 12.1 (Core null form estimates). The following hold: 

(131) || PkM^k^S) Ik < 2 fcl 2^- max ^^}) 2 ^| fc2 - fe 3||| (2) w (3) || fli ^ 

(132) || (/ - KMtfM?) Ik < 2 fcl || 0« || 5 x|| < } ||^, h < h - C , 

(133) || n^^uh^S) WlHV) ^ 2fcl ll 4$ llz^H 4? ||*, fci < - C , 

(134) IKi-^JP^^,^)!!^ 

< 2 fc i2 5 ( fcl - max ^ 2 ' fc3 »2" <5|fc2 ~ fc3| || 0^ || s i|| 0^ || 5 i , 

(135) \\n k M^,^)\U^ 

< 2 k ^ k ^U^\\ sl \\<f ) ^\\ sl , h > max{fc 2l k 3 } - C . 

In addition one has the following quadrilinear form bounds, which hold under the condition 

(136) K^n^ ■ d a ^),n k (d a ^ ■ «>| 

< 2^ ^})||^)|| 51 ||^)|| 51 ||^)|U 1 ||^ 4 || JV ,, 

(137) KCDA)- 1 ^^ • d^),d t H k {d t ^ ■ M)\ 

< 2^*>)||^||5 1 ||^ ) ||^||^||5i||^||^, 

(138) |<(DA)- 1 V :E ^(0l 1 i ) • V^J.M^-^))! 

< ^^Mll^ll^ll^ ||5X||^||^||^ 11^. 

The second category of bounds is covered by: 
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Theorem 12.2 (Additional core product estimates). The following hold: 

(139) \\(I-H* kl )(Q <k2 - C A kl d t K)\\N < 2l fcl ||A fcl || i2(L2) ||0 fe2 || 5 i, h<k 2 -C , 

(140) || W kl (A kl d t (t) k2 ) \\l^(l 2 ) ^ II \\z«u\\ <t>k 2 Ws 1 , h<k 2 -C, 

(141) || (/ - n^P^dt^S) \\az* < 2 5(fel " fc2) || 0g 3 Ml C ||* , h<:k 2 -C , 

We conclude this Section with the application of Theorems H2.ll and ll2.2l to the estimates 
listed above. This is for the most part straight forward and left to the reader. 

Proof that TheoremUHl implies estimates (|3"I]), (|32|), (153]), (|53j), (155]). and (|5D|). For ([31]) we 
use (II5Tj> . for (J32J) use (HMD, for ([53]) use (JT3T]) with k 2 >k x + 0(1), for (JMD use ([132]), for 
(ESD use (1153)1 . and for QSDD use JSSJ) - (fI58|» . □ 

Proo/ £/ta£ Theorem\7£M implies estimates (]42]) ; ([58]) , and ([59]). For g2]) use (1T4T]) . for §E> 
use (TT39]) . and for ([59]) use (IT40D . 

□ 

For the remainder of this Section we prove Theorems 112.11 and 112.21 We begin with a 
simple calculation that will be used a number of times: 

Lemma 12.3 (Square summed L 2 (L 6 ) L 2 (L°°) estimate). Letj-C ^ k'+l' ^ ±(j+fc)+0 
with V < C and k' ^ k + C . Then one has the following uniform estimate: 

(142) ( \\Pc k ,(i>)Q< 3 <Pk\\h {L ^ < 2l fc '2^'2ll0 fc || 5fc[L2(L6)] , 

Cy{l') 

where Cy (I') is a finitely overlapping collection of radially oriented rectangles of dimensions 
(2 k +l ) 3 x 2 fc . Here S k2 [L 2 (L 6 )] refers to the L 2 (L 6 ) portion (including square sums) of the 
norm from lines and ([S]) • 

Proof of estimate (11421) . This follows immediately from the L 2 (L 6 ) estimate on line ([H]) and 
Bernstein's inequality in the form Pc k ,(i>)L 6 C 2^ k '2^ 1 ' L°°. □ 

12.1. Proof of the null form estimates. We begin with estimates ( 1 1 3 II) — ( fT33l) . These 
can be boiled down to an even more atomic form which is the following: 

Lemma 12.4 (Core modulation estimates). The following estimate holds uniformly in the 
indices ji, ki, where j 2 ,j3 = ji + 0(1): 

(143) | {Q h </$ , Af (Q<j 2 4> ( k ] , Q <j3 (f> { ^ ) > | 

< a-^-fc^-^.-fcai^^a}^!! a) u u (2 } || n (3) hoo{L2) ^ 

In addition, when j > k m i n + C one has the improved bound: 

(144) | {Q^M^fl <2"«y- fc -)2 2fc -2*»- \\<j>% || o,, ||< } || 5fc2 \\n* • 

Proof of estimate ( I143j) . There are three cases depending on the relative separation of spatial 
frequencies. Each of these is further split into low and high modulation subcases. 
Case la:(&2 = k$ + 0(1) and j\ < k\) Here the frequency angle separation between the 
second two factors is Z(0 2 ,</> 3 ) < 2~ k2 2^ kl+ ^' := 2 l . The null form M saves one power 
of this angle. On the other hand, the angle of interaction with the output is Z^ 1 , <ft 2 ) < 
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22O1 := 2 1 ' . Breaking up the high frequency factors in the null form with respect to the 
2 kl radial scale and the 2 l sector scale it remains to estimate the expression Qj ± F in L 2 where 

F:=P kl ^(P^OQ^S^-C^^Q^) 

C H {V) 

Disposing of M we start with the fixed time estimate 

c kl (i') 

< 2 2 k3+kl+n I £ \\P Ckl{n Q <h cf>W(t)\\lA I £ \\P-e kl (i')Q<^S(t)\\^ 
\c kl (i') J \c kl (i') 

where at the last step we have used orthogonality in frequency and then the L°°L 2 bound- 
edness of Q<j 3 - Hence integrating in time we arrive at 

(145) ||F||£ 2 < 2 2fc ^i0g)||2 ooL2 \\Pc kl{n Q<^(t)\\l^ 

c kl (i) 

Then applying Lemma [12.31 to put 0? in L 2 (L°°), we see that: 

\\QnF\W <2§^<'2§ fc ^+^>|| $ |k 2[W)] || 0® 

= 2 lh 2 \^) 2 W^)2 k -2 k ^\\ <f>® \\s k2[ mm\\ <f>% 
concluding the proof in this case. 

Case lb:(A;2 = k$ + 0(1) and ji ^ k\) Here integrating by parts we get a factor of 2 kl+k2 
from the null form. Then we localize 0j^ and 0j^ with respect to 2 kl sized frequency cubes, 

but without any modulation localization. For 0^ we use the bound (jHJ), and for 0j^ we use 
the L°°L 2 norm. Then the same computation as above for 



F = ^{Pc kl Q<n<P k ^P-c kl Q<rJ2) 
yields 

||F|| L2 <2^+l fc i0i 2 2 )|| Sfc2 ||0£)|| iooi2 

which suffices. 

Case 2a:(fci = k 2 + 0(l) and ji < k 3 ) This case is mostly analogous to Case la. Integrating 
by parts one may place the null form between 0^/ and 0^ . Then the angular decompositions 
between the inputs of these two terms, the angle of the output (frequency 2 k ' A now) is the 
same as Case 1 but with k± and k% transposed. With 

F-.= p k:i £ ^f{Pc ks (nQn<l> k 1 ^P-c k3 (nQ<^) 
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the relevant computation becomes: 



C k Jl') 



1 1 

2 / \ 2 



c* 3 0') / \c fc3 (i') 



< 2 ii 12 i(i 1 ^3) 2 i(^- fc2 ) 2fcl2fc22fc3|| g . i0 (l) ||l2(l2)|| ^(2) || 5fc2[L2(i6)]; 

which suffices. 

Case 2b:(A;i = &2 + O(l) and j\ > k 3 ) The same argument as in Case 2a applies, with the 
only difference that cp^ and 0^ are now decomposed on the frequency scale 2 fc3 , and there 
is no modulation cutoff. Thus the bound has to be used for the latter. 
Case 3a:(fci = k 3 + 0(1) and j < /c 2 ) The computation here is essentially the same as 
above, but with slightly different numerology. The main difference is now that Z(</> 2 , <f) 3 ) < 
22^ 1_fc 2) := 2 l , and we sum over sectors of size ~ 2 k2 x (2 fc2+i ) 3 . Then a computation similar 
to Case 1 gives us: 

LHStfUSD < 2-^2i l 2f*»2*»|| $ \\ ^ \\ <fij> lk 2 [ iW] || J> { £ , 

_ 9i(ji-fc 2 )92fc29fc3|| A 1 ) II || J 2 ) || „ „ || A?) || 

- > 2 2 || (p ki || 0i i|| (p k2 \\s k2 {Li(L0)}\\ fc3 • 

Case 3b:(fci = &3 + 0(1) and j > k^) Here we directly use the easy product estimate: 
LHS(H13} < 2-^2^2 fc 3||0«|| 0ii || ^ || Sfca || ^ |U- (La) . 

□ 

Proof of estimate (11441) . There are two cases depending on the relation of the modulations 

of 0^ and </>® to 2K 

Case l:(One of or 0^ /ias modulation comparable to 2 J ) This situation is symmetric. 
It suffices to show the single estimate: 

llPkMQ^-c^,^)^^) < 2^2^^^\\^ || nill^lU-^) , 

which follows immediately from putting a derivative of M on the lowest frequency and 
L 2 — > L°° Bernstein's inequality for the lowest frequency as well. 

Case 2:(Both </>£f and (j)£ have modulation <C 2 J ) This can happen only when k\ < /c 2 — C, 

in which case it is a (++) or ( ) type High x High =^> Low interaction between <jy£' and 

(jyj^ ■ Thus, j = k 2 + 0(1) and it suffices to prove: 

llPkMQ^-c^^Q^-c^h^) < 2^2^11^11^ || *2 } IU-C£») • 

This follows at once by putting a derivative on the low frequency output, and breaking the 
product into antipodal blocks of scale 2 fcl while using the special L 2 (L°°) square sum 
bound © for the first factor. Note that the multiplier Q K j-cPc k is disposable thanks to 

j > h + C. □ 
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Proof of estimate (I13ip . This follows from the next two bounds: 

(146) || (/ - n kl - Q >kl+c )PkM<P { ^ lUw S LHSdmD , 

(147) || (Q>k 1+ c + n kl )p k M^^ { S) ILo,-i < LHSdm]) . 

Case l:(Estimate (11461) ) The restriction induced by (I—'Hk 1 — Q>ki+c) means that one of the 
two input factors always has the leading modulation. First permute notation so the output 
frequency is k% and the two inputs are k±, hi. Then by duality (note that M is skew-adjoint 
as a form) and symmetry of the estimate, it suffices to show bounds of the form: 

\(Q j( f)^ ,M(Q <j+ o(i)<f>k^Q<j+o(i)^k 3 ))\ 

< 2 -<5| J -A ;2 | 2 A : 3 2 <5CA : 3- m — {^x,A3 2 » 2 -,5|A :i _ fc2 | 2fcl2fc2 1 1 ^£L) 1 1 ^ ^ 1 1 <^£^ 1 1 ^ fc2 || ^ 3 IU~(2*) , 

which follows directly from (j!43p . 

Case 2a: (Estimate (I147P for low modulations) Freezing the output modulation it suffices to 
show: 

II Q 3 S kl N{Q <3 -c4> { ^Q< 3 -c<i>f s ) II n,_ 4 

< 2~ 5 l^~ fe2 l2 fel 2^ fcl ~ max ^ fc2 ' fc3 ^2 _5 ' fe2 ~ fc3 '2 fc2 2 fc3 \\d>^ II o \\d>^ II <? 
which again is a direct consequence of (11431) . 

Case 2h:(Estimate (I147P for high modulations) In this case we use: 
WQ^U^t^S) II 0,-1 < 2~^^- fc -l2 fc -2^^~— ^- fe ->>2-^l^-^l2^2 fc -||0g ) |U fc2 ll^ll^ , 

follows immediately from (I144p in the case j > k\ + C . □ 



Proof of estimate (I132p . There are two main cases: 

Case l:(</4^ with highest modulation) Due to the restriction oiTi*, this can only occur when 
the second factor is Q>k x +c ( Pkl '■ Then we use (I144p which directly implies for this case: 

\Wm { £J£)\U < (2^ + 2^1)2^11^1^11^11^ • 

Case 2:(Output or (j^Q have the leading modulation) In this case we end up needing bounds 
of the form: 

|| QMQ<^o { i)^},Q< j+ o { i)<l> { £) IL o,- 4 < 2-^l2^2 fe2 ||0W|| 5fci ||^ ) || 5fc2 . 
IIQ^ow^g^ow^^Q^)!!^^) < 2^^l2 2fc 1 2 fc2 ||0i 1 i ) || 5fci ||0g ) || 5fc2 , 
which are both immediate from (11431) . □ 



Proof of estimate (I133P . This is a direct computation using angular decompositions. At a 
fixed modulation 2 J the angle of interaction is Z(0 1 ,0 2 ) < 2^~ kl ^ := 2', and disposing of 
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the null form we have: 

< 2 l 2 kl 2 k2 II PiQki+2i<i>k 1 ■ II PiQ<j-c4>k 2 IU°°(i 2 ) • 

OJ 

Using Cauchy-Schwarz and: 

II P\Q<3-c4k] IIl°°(L2) < ^ II P 2^-k 2 )Q<3-C<f>k2 lll°°(i 2 ) ' 

we have: 

II Q^-cM^l Q <3 -c4> { £) < || 0« || <w || \\ sl2 , 

which suffices. □ 



Proof of estimate (I134p . This follows immediately from (j!4ip and ( fTTJj) . Notice that the lack 
of an L l (L 2 ) estimate for Q>k 1 +cPk^{ ( j>k} 'j^fe) i s irrelevant because definition of Z hyp 
limits modulations to <5<fc 1 +c-Pfci- D 

Proof of estimate (11351) . By symmetry we may assume /C2 < £3 + 0(1), in which case k\ = 
fa + O(l). There are two subcases depending on the size of the output modulation 2 J : 
Case l:(/c2 > j—C) In this case the angle of interaction between the two inputs is ^(</> 2 , 3 ) < 
2|0'-fc2) ._ 2 1 ' '. On the other hand, the output sector localization of 2' := 2^~ fel ) « 2^ j ~ k,i) 
is passed to the high frequency factor, so using Lemma [12.31 we have: 

2ll PtQjSkMQKj-C^ktQ^-C^S) \\l H L~) 

< 2'2 2 ''2^2 2fe 3 £ £ || P^Q^® ||i w || P^-^i? Hi,^, , 

< 9 2« 9 3i' 9 5fc295fc3|| J. ( 2 ) ||2 || j,(3) i|2 

^ 2 2 2 2 || </) fc2 || Sfc2 [ £ 2( L 6)]|| ^ fe3 HSfcg [£»(£«)] > 

Case 2:(/c 2 < j — C) Here the calculation is essentially the same as above, except that V = 0. 
This again gives: 

|| QjSuMQ^-c^hQ^-c^) < (2 J+fcl )-2 fcl 2^- fc 3)||0^|| 5 ,J|0g)||^ . 

□ 



Proof of estimate (11361) . Freezing the output modulation, we will show: 

(148) | (□- 1 Q J P fe (Q <J „ c 4 1 i ) ■ daQ^c^l QM^Q^-c^S ■ Q<j-cM)\ 

< 2 iW-*) 2 §(*-^})|| 0« || sl || 0g) || 51 || 0g) || sl || ^ \\ N * . 

Here we are in the configuration k\ = /c2 + 0(1), = k^ + 0(1), and k < ki — C . Thus, 
the left and right products are summed over diametrically opposite angular sectors of size 
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2 k x (2 fc+i ) 3 , where 2 l := 2^~ k \ On the other hand, the null form between the second and 
third terms gains us Z(0 2 , 3 ) 2 , where this angle mod n cannot exceed 2 l+c . Therefore we 
group the product of the two diagonal sums into dyadic values of Z(0 2 , 3 ) mod 7r and break 
into three cases: 

Case 1:(Z(0 2 , 3 ) mod 7r < 2 l 2 k ~ k2 ) A little care is needed to use orthogonality in space. To 
gain this, at first keep the second diagonal sum under the time integral as follows: 

LHS([148J) |z(02^3) mod 7r<2 ( 2 fe - fc 2 

< 2 -k 22k:i || P Ck{l)Q<j _ c ^ || i2(ioo) || P^Q^-a^S !Uw 

x sup || Pc^oQ^j-c^C*) IU||| P-c k {i)Q<i~c^{t) \\li ■ 

The inner sum of the second factor on the RHS can easily be reconstructed after Cauchy- 
Schwarz by spatial orthogonality. On the other hand, for the two L 2 (L°°) norms we use 
Lemma [12. 3 1 This gives us: 

LHS ([HE]) I z( 02 ^3) mod 7r < 2 ; 2 fc-fe2 

< 2 |^ 2 | fcl2 -i fc22 , 3 || (i) || Sfci[i2(L6)] || g) || 5fc2[L2(i6)] || 4^ \\ Sk3 1| ^ \\ N , , 

< 2^-*) 2 S(*-^*»|| 0« IMI $ \\s4 ^ Ml ^ 4 Ik* , 

which is even better than (I148p . 

Case 2:(Z(0 2 , </> 3 ) mod 7r < 2 l 2 k ~ k3 ) This is essentially the same as Case 1 above, but since 
the angular gain is in frequency 2 fcs we put 0^ in a square summed L 2 (L°°) via Lemma Tl 2. 31 
and (f)^ in L°°(L 2 ) instead. This gives: 

LHS ([USD |/(02 ^3) mod7r < 2 i 2 fc-fc3 

< 2 l fc 2'2^2f^2-^3|| W ||^ (L2) || g) || 5fc2[L2(L6)] || 0g) || Sfc3[i2(i6)] || \\ N , , 

< 2 i(^,) 2 !(^n { M)|| « || 51 || g) || sl || g) W ^ ||^ j 

which suffices. Note that a similar square summation for fixed time procedure as in Case 1 
was used here. 

Case 3 : (2'2 fe-nun { fe2,fe3 } -C Z(0 2 , 3 ) mod n < 2 l ) In this case there is a definite angle between 
spatial frequencies in the blocks Ck{l),C' k {l) for fixed dyadic 2 l = Z(0 2 ,0 3 ). Thus, we can 
use one power of Z(</> 2 , <f) 3 ) to put the product of <f>^ and 0^ in L 2 (L 2 ) using null frames. A 
little more care is needed here to gain spatial orthogonality for ip^- Thus, we first fix time 
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and compute: 

LHS Q148J) | z ^ 2 ^3) mod ^^2'' ,t=const 

< 2 -2( k+ i) 2k22 k 32 2i' j2 || PcmO^-c^W • P±c> k (i)Q<^ { £(t) \\ Ll 

C k (l),C' k (l): 
Z(C fe (0,±C fc (/)')~2 ; ' 

X || Pc k {t)Q<j-C<l>k?{t) |U«|| P^C' k {l)Q<j-C^kM lUi > 

< 2 -2( fe+ o 2fe22fc32 2r|| g< ._ cV , fc4 ( t) | Ur ( £ || P^Q^V® HI.) 1 

x ( E II p cmQ<i-o^^) ■ P*wQ<i-c4>$V) llii) 1 • 

Ci(0: 
^(C fe (0,±C fc (Z)')~2 i ' 

Integrating and using Cauchy-Schwarz in time, and then using the special microlocalized 
L 2 (L°°) block norm from line flS} for 0^, we get: 

LHS(H48])| z( , 2 ^ )mod _ 2 , < 2- fc 2- 2 ¥'2^2^2 fc 3J 23 (r)||0«|| 5fci ||^ 4 ||^ . 

where: 

/ 23 (/') 2 = E 2 2 m P Cfc (oQ<i-c0i 2 2 } • P*wQ<j-o<l>$ \\ 2 »M ■ 
c k (i),c' k (iy. 

Z(C fc (0,±C fc (0')~2 ! ' 

Since we already spent 2 fc3 we put the second factor in L^(L 2 ± ) and use L 2 (L°^) for the 
first. This gives us: 

and so: 

LHS(HlHl)| z(02i ,3 )mod _ 2 , < 32^2*^-— 0£ 5 M^ll 0g> H^xll H^ll V'** IIJV • 
Summing this over all I' < I + C gives (j!48p for this case. □ 

Proof of estimate ( I137P . Freezing the output modulation, we'll show: 

< S2*<^— *>22*C*— *->|| H^ll ^&«> H^ll . 

We will use an L°°(L 2 ) estimate for both the third and fourth factors. Thus, via Holder's 
inequality, a simple weight calculation, and the Leibniz rule we have reduced matters to: 

(149) || QjP^Q^c^ ■ tPQ^c^S) \\l H l~) < V+*2ft- k W-V\\ 0« M| <t> { £ \W , 

(150) || QMQ<i-c^ ■ □Q<i-c^ 2 2 ) ) \\lHl~) < 2^2§t>- fe >2^)|| 0« || sl || <f><£ ||* . 

To prove (I149p . note that the null form gains two powers of the angle Z(0 X ,0 2 ) < 
2 2( k+ 'j)2~~ kl ■= 2 l . Therefore, breaking the product into a sum over antipodal radially di- 
rected blocks of dimension 2 k x (2 k+l ) 3 , where I' = ~(j — k), and using (11421) for both factors 

45 



we have: 

LHSdHi < 2 a 2''2^2^2^||^ 1 i )|| Sfci[L2(L6)] ||^) \\ aka[ » m , 

< 2^2^^)2l^)2^2 fc l< ) || 5fci ||< ) || 5fc2 . 

To prove (11501) we use the same calculations as above, except for the second factor we 
trade flUgJ) for: 

(151) ( £ II Pw^-o^ lli 2 ^)) 1 < 2 2fc 2§<'2^2 fc l ^ || j , 

cm 

Q 1 j 

which is an immediate consequence or Bernstein's inequality and OQ < j_ c P k2 X ^ 2 C 2^2 k2 L 2 {L 2 \ 
This gives: 

LHS(H50D < 2 j+2k 2^ j ~ k h^ k - kl) 2 kl 2 k2 \\(f)ff \\s kl \\(f>^ lk 2 ■ 

□ 

Proof of estimate (I138p . Freezing the output modulation, our goal here is to show: 

< 2i«-*) 2 *c*-^»|| C ii^n M| C |M| ^ 4 Ik • 

Expanding the null form and computing the weights, it suffices to have: 

(152) || QMQ^-c^} ■ VxQKj-ctfg) < 2i*2iy-«2i(*-^|| || 5 i|| < } || fl i , 
(153) 

II QMdaQ^-c^ ■ d^^cM < 2^2l fc 2i^- fc )2^ fc - fe 3)|| ^jj || sl || ^ \\ N , , 

(154) || QjPkiDQ^ct® ■ Q^-c^) \\v>w < y2i*2^'-*)|| ^ || s i|| ^ 4 Ik* • 

The proof of these three estimates follows from the same angular decompositions and antipo- 
dal block sums as in the last proof. Estimate ( 1152[) follows by using (11421) for the first factor 
and L°°(L 2 ) for the second. Estimate (I153j) follows similarly once once the null form is taken 
into account. Note that here the antipodal block sum for ipk 4 needs to be reconstructed for 
fixed time using orthogonality. Finally, estimate (11541) is the same as (I153f) but uses (11511) 

(3) 

instead for (J) h \ Further details are left to the reader. □ 

12.2. Proof of the additional product estimates. These estimates use the same kind 
of modulation/angular-sum decompositions as in previous proofs, so we leave more of the 
details to the reader. 

Proof of estimate H139[) . There are three main cases: 

Case l:(^4fei with highest modulation) Due to the restriction of H*, this can only occur when 

the second factor is Qki+c< <k 2 -c^k 1 - Then there are two subcases: 

Case l&:(Contribution of dtQ > k 1 ~c4>k 2 ) Here we use a product of the two bounds: 

II Qk 1 +c< <k 2 ~c^k l ||l 2 (l°°) ^ 2 2fcl || A kl ||l 2 (l 2 ) , II d t Q >kl -c4>k2 ||l 2 (l 2 ) ^ 2~i kl \\ d t <f>k a \\ ,\ ■ 

Case lh:(Contribution of dtQ<k 1 -c4 > k 2 > an d output modulation > 2 kl ) Note that this is the 
remaining case because an output modulation of <C 2 kl is impossible due to the restrictions on 
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Qk 1 +c<-<k 2 -c^-k 1 and d t Q <kl -c4 ) k2- Here we use a product of the previous L 2 (L°°) estimate 
for the first and L°°(L 2 ) for the latter. 

Case 2:(0utput or <p k2 have the leading modulation and j < k\ + C) In this case we can 
reduce things to the bounds: 



\\Q j (Q<j+o(i)A hl d t Q <j +o (1) <l> ha )\\«,-l < 2^- fel )2t fcl 2 fc2 || A kl IU^jH^ 

x i 

Q <J+0{ i)(Q< J+ o(i)^d t Q 3 ^) || il(L2) < 2^hl k ^\\ A kl || L2(L2) || <p k2 1| 



Both of these bounds follow from the usual angular sum decompositions and the Bernstein 
type estimate: 

( E II P c kl (^))Q<j-oA kl < 2^2i^|| A kl || L2(L2) . 

C H (\(j-ki)) 

Case 3:(Output or <p k2 have the leading modulation and j ^ k\ + C) This is analogous to 
the last case except there are no angular decompositions and instead we simply use: 

II Q<j-cA kl ||l 2 (l°°) < 2 2fcl || A kl \\l 2 (l 2 ) , 
instead of the previous square sum bound. □ 

Proof of estimate (1140}) . This follows at one from the definition of Z el1 , which allows one to 
sum over the modulation localizations enforced by ^H kk - D 

Proof of estimate ( I14ip . Freezing the output modulation of the product, and getting rid of 
the contribution of H kl , it suffices to show: 

|| QjP^Ta |Ui(L») < 2^2^^)2i( fcl - fe )||0g ) || 51 ||0g ) || 51 , 
for either one of: 

Both cases are essentially the same due to the matching fc 2 = k 3 + 0(1). After breaking the 
sum into angular sectors one uses: 



;(2) 



; fc 2 ll„o,i , 
C fcl (i') 

for the high modulation factor and (1142 j) for the low modulation factor, both with V = 
|0 — ki)- Note that can be further broken into high and low modulation pieces, and 

the product of two high modulations is even more favorable. □ 

Appendix: A multilinear null form for MKG-CG 
Schematically we have: 

DAi = -P^(0V0) , AA = -3(05^) , 

where V is the Hodge projection. In 3D one can write this operator conveniently in terms 
of the vector product Vx. In higher dimensions it is easier to use the simple formula: 

V = I- VA-V , 

where the first V is grad and the second is div. 
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Now investigate the expression: 

N{A,4>) = A a d a (f) . 

Directly plugging things in we have: 

QiQi 

Af(A, <p) = A~ l %<pd t <p) ■ d t <p - a-^d^) ■ + —Sttfdit) ■ d j( j) . 

To uncover the null structure, we add and subtract the expression □ _1 S(0<9t0)<9t0. This 
gives us: 

Af(A,<j>) = -Q 1 (A,<f>)+M 2 (A, ( t>) , 

where: 

Af 2 (A,<p) = A^3(^^) • 9^ - ^) • 9^ + • 3> . 

To continue the computation we use A -1 — CT 1 = — A □ so the second term simplifies 
to: 

tf 2 (A, 0) = -A^D-^^dtcf)) ■ dt<f> + A^a^d'd^^d^) ■ drf . 
This can be further reduced to the sum of Qo null structures by adding and subtracting (for 
instance) d l dt^s(4>di<j))dt4>, which gives: 

= Q 2 (A0) + Q 3 (A0) , 

where: 

Q 2 (A0) = A- 1 ^- 1 ^^^)-^, 

Q 3 (A,cf>) = A-in-^&SSi^ .&*<!> . 

There is a conceptual way to visualize this by duality. The expression A a d a cj) is a sum 
of three Q null structures: The interaction is between the first two factors, the second two 
factors, or mixed between 2 and 3 . The latter is the main one which involves null frames. 
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